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SUMMARY 


A  method  for  predicting  the  static,  longitudinal  aerodynamic 
characteristics  of  typical  missile  configurations  at  zero  roll  eingle 
(i.e.,  in  a  plus  configuration)  has  been  developed  smd  programmed  for 
use  on  the  IBM  7090  digital  computer.  It  can  be  applied  throughout 
the  subsonic,  transonic,  and  supersonic  speed  regimes  to  slender  bodies 
of  revolution  or  to  nose-cylinder  body  combinations  with  low  aspect- 
ratio  lifting  .surfaces.  The  aerodynamic  characteristics  can  be  com¬ 
puted  for  missile  configurations  operating  at  angles  of  attack  up  to 
180  degrees.  The  effect  of  control  surface  deflections  for  all  modes 
of  aerodynamic  control  are  token  into  account  by  this  method.  The 
method  is  based  on  well-known  linear,  nonlinear  crossflow  and  slender 
body  theories  with  empirical  modifications  to  provide  the  high  angle 
of  attack  capability.  Comparisons  of  the  theory  with  experimental  data 
are  presented  to  demonstrate  the  accuracy  of  the  method. 
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NOMENCLATURE 


AR 

b 


d 


f 


i 


exposed  aspect  ratio 

semispan  of  an  aerodynamic  surface  including  the  body 
radius,  feet 

total  drag  coefficient 

base  drag  coefficient 

crossflow  drag  coefficient 

friction  drag  coefficient 

induced  drag  coefficient 

total  zero-lift  drag  coefficient 

pressure  drag  coefficient 

wave  drag  coefficient 

incompressible  skin- friction  coefficient 
compressible  skin- friction  coefficient 
drag  of  a  flat  plate  normal  to  the  flow 
total  lift  coefficient 
lift  curve  slope,  per  radian 

total  longitudinal  pitching  moment  coefficient 
root-chord  of  an  aerodynamic  surface,  feet 
tip-chord  of  an  aerodynamic  surface,  feet 
diameter  of  the  body  at  any  station,  feet 
base  diameter  of  the  body,  feet 

diameter  of  the  nose  at  the  nose-body  Juncture,  feet 

spanwise  location  of  the  vortex  which  emanates  from  the 
forward  surface,  feet 

height  of  the  trailing  vortex  above  the  body  centerline 
at  the  aft  surface  center  of  pressure,  feet 

downwash  interference  constant 
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NOMENCLATURE 

(continued) 


apparent  mass  factor 

linear  lift  interference  factor  due  to  angle  of  attack 

linear  lift  interference  factor  due  to  control  surface 
deflection 

total  length  of  the  body,  feet 
length  of  the  nose,  feet 

arbitrary  reference  length,  usually  the  maximum  body 
diameter,  feet 

distance  from  the  tip  of  the  nose  to  the  intersection  of 
the  tail  leading  edge  with  the  body,  feet 

distance  from  the  tip  of  the  nose  to  the  intersection  of 
the  wing  leading  edge  with  the  body,  feet 

free-stream  Mach  number 

cotangent  of  the  leading  edge  sweep  angle 

radius  of  the  body  at  any  station,  feet 

Reynolds  number 

base  area  of  the  body,  (feet)^ 

exposed  planform  area  of  one  pair  of  forward  lifting  sur¬ 
faces,  (feet)^ 

body  cross-sectional  area  at  the  nose,  body  juncture,  (feet) 
planform  area  of  the  body,  (feet)^ 
surface  area  of  the  body,  (feet)^ 

exposed  planform  area  of  one  pair  of  tail  surfaces,  (feet)^ 

planform  area  of  one  pair  of  tail  surfaces  as  obtained  by 
extending  the  leading  and  trailing  edges  to  the  center- 
line  of  the  body,  (feet)^.  See  Figure  2. 


NOMENCLATURE 

(continued) 

exposed  planform  area  of  one  pair  of  wings,  (feet)^ 

planform  area  of  one  pair  of  wings  as  obtained  by 

extending  the  leading  and  trailing  edges  to  the  center- 
line  of  the  body,  (feet)^.  See  Figure  2. 

volume  of  the  body,  (feC^)^ 

distance  from  the  nose  to  the  missile  center  of  gravity, 

( feet )  ^ 

distance  to  the  surface  center  of  pressure  as  measured 
from  the  intersection  of  the  leading  edge  of  the 
aerodynamic  surface  with  the  body,  feet 

distance  from  the  nose  to  the  center  of  pressure  location, 
feet 

distance  from  the  intersection  of  the  panel  leading  edge 
and  the  body  to  the  hinge  lipe,  feet 

distance  from  the  nose  to  the  centroid  of  the  body  plan- 
.  form  area,  feet 

missile  angle  of  attack,  degrees 

compressibility  factor,  |m^-1 

control  surface  deflection,  degrees  (See  Figure  I  for  sign 
conventions ) 

component  of  the  induced  drag  coefficient 

increment  of  wave  drag  for  the  transonic  speed  regime 


ratio  of  the  drag  coefficient  of  a  circular  cylinder  of 
finite  length  to  that  of  infinite  length 

conical  nose  semi-vortex  angle,  degrees 

lifting  surface  taper  ratio,  C  /C 

leading  edge  sweep  angle,  degrees 

sweep  angle  of  the  quarter  chord  line,  degrees 


NOMENCLATURE 

(continued) 


SUBSCRIPTS 

A  aft  lifting  surface,  alone 

B  body  alone 

BT  body  in  the  presence  of  the  tail 

BT  -  a  body  in  the  presence  of  the  tail  due  to  angle  of  attack 

BT  -  6  body  in  the  presence  of  the  tail  due. to  control  surface 

deflection 

BW  body  in  the  presence  of  the  wing 

BW  -  a  body  in  the  presence  of  the  wing  due  to  angle  of  attack 

F  forward  surface  alone 

FB  forward  surface  in  the  presence  of  the  body 

N  nose 

T  tail  alone 

T  -  a  tail  alone  due  to  angle  of  attack 

T  -  6  tail  alone  due  to  control  surface  deflection 

TB  tail  in  the  presence  of  the  body 

TB  -  a  tail  in  the  presence  of  the  body  due  to  angle  of  attack 

TB  -  6  tail  in  the  presence  of  the  body  due  to  control  surface 

deflection 

TV  ,  tail,  nonlinear  component 

W  wing  alone 

WB  wing  in  the  presence  of  the  body 

I 

Vffi  -  a  wing  in  the  presence  of  the  body  due  to  angle  of  attack 

WV  wing,  nonlinear  component 

The  control  surface  is  defined  as  the  tail  regardless  of  the  mode  of  con¬ 
trol;  the  fixed  surface  is  defined  as  the  wing  (see  Figure  l). 
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INTRODUCTION 


Increasing  maneuverability  requirements  of  missiles  indicated  a 
need  for  predicting  the  aerodynamic  characteristics,  including  lift, 
drag,  and  pitching  moment,  of  missile  configurations  to  angles  of 
attack  of  90  degrees  and  higher-  A  study  shoved  that  existing  methods 
for  computing  these  aerodynamic  characteristics  are  based  on  a  number 
of  different  theories  all  of  which  are  applicable  only  to  small  angles 
of  attack.  To  fulfill  the  high  angle  of  attack  requirements,  a  method 
for  determining  the  aerodynamic  characteristics  of  low  aspect-ratio 
configurations  at  zero  roll  angles  operating  at  angles  of  attack  up  to 
180  degrees  has  been  developed.  The  method  is  applicable  throughout 
the  subsonic,  transonic,  and  supersonic  speed  regimes  up  to  g'AR  ^  10-0, 
and  accounts  for  control  surface  deflections , 

The  method  is  composed  of  well-known  linear,  nonlinear  crossflow, 
and  slender  body  theories  which  have  been  modified  to  provide  the 
required  high  angle  of  attack  capability.  These  theories  can  be 
applied  to  slender  bodies  of  revolution  or  nose-cylinder  bodies  with 
canard,  wing,  or  tail  controls  (Figure  l). 

This  report  describes  the  methods  developed  and  the  computer 
program  which  has  been  written  for  use  on  the  IBM  JO9O  digital  computer. 
The  description  of  the  method  is  divided  into  three  parts;  lift,  drag, 
and  pitching  moment.  For  the  sake  of  clarity,  the  description  of  the 
method  is  kept  to  a  minimum,  without  lengthy  justification  and 
description  of  the  techniques  employed*  The  reader  is  referred  to  the 
references  for  detailed  descriptions  of  the  various  theories.  The 
description  of  the  computer  program  consists  of  a  brief  discussion  of 
the  main  program  and  subroutines  and  complete  instructions  required  for 
use  of  the  program.  Comparisons  of  theoretical  results  with  experimental 
data  are  presented  for  angles  of  attack  up  to  90  degrees  over  the  entire 
speed  range  to  demonstrate  the  accuracy  of  the  theories.  Some  data  for 
a  missile  configuration  at  I80  degrees  angle  of  attack  is  available  and 
is  compared  with  the  theoretical  results, 


LIFT  CHARACTERISTICS 


The  total  lift  on  the  missile  is  the  sum  of  the  body  lift,  the  lift 
due  to  the  aerodynamic  surfaces,  and  the  interference  lift  between  the 
forward  and  aft  surfaces.  The  lift  on  the  body  and  aerodynamic  surfaces 
is  composed  of  two  components;  linear  lift  including  the  effects  of  tlw 
body-lifting  surface  interaction  and  nonlinear  crossflow  lift.  In 
general,  the  crossflow  lift  component  is  caused  by  flow  separation  whf® 
occurs  at  angle  of  attack,  while  the  interference  component  is  the  lift- 
loss  on  the  aft  lifting  surface  due  to  downwash  from  the  forward  surface 
(Reference  l). 

Allen,  References  2  and  3,  developed  a  method  for  predicting  the 
total  lift  on  bodies  of  revolution  at  angles "of  attack.  This  method 
includes  the  linear  or  potential  flow  component  and  two  nonlinear 
components:  the  viscous  crossflow  force  and  the  viscous  axial  force. 
Because  the  contribution  of  the  social  force  component  to  the  body  lift  * 
is  small,  it  is  usually  neglected.  Allen’s  expression  for  the  body 
lift  is 


where  the  first  term  is  the  linear  contribution  and  the  second  term  is.  the 
nonlinear  contribution.  The  appareht  mass  factor,  k^  "  k^,  and  the  drag 
ratio,  n,  can  be-  obtained  from  Figure  3,  while  the  crossflow  drag 
coefficient,  is  obtained  from  Figure  1+.  Comparisons  of  theory  with 

experimental  data  for  numerous  bodies  of  revolution  over  a  wide  range 
of  Mach  numbers  and  angles  of  attack  are  presented  in  Reference  3.  It 
should  be  noted  that  although  this  expression  for  the  lift  is  independent 
of  the  nose  shape,  good  agreement  with  experiment  is  indicated  in 
Reference  U  for  a  body  with  an  unu^^al  shape. 

The  linear  lift  characteristics  low  aspect-ratio  lifting  surfaces 
whose  cross-sections  are  thin  and  symmetrical  are  generally  a  fvinction  of 
speed,  planform  area,  and  aspect-ratio.  When  the  diameter  of  the 
missile  body  is  of  the  same  order  of  magnitude  as  the  span  of  the 
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lifting  surfaces ,  the  effects  of  hody-wing  and  "body-tail  interactions  are 
significant.  Hence,  the  linear  lift  of  the  aerodynamic  surfaces  is  composed 
of  two  components:  the  lift  on  the  surface  in  the  presence  of  the  body, 
and  the  added  lift  on  the  body  due  to  the  presence  of  a  surface.  Most  low 
aspect-ratio  missile  configurations  exhibit  a  nonlinear  dependence  of  lift 
on  angle  of  attack,  especially  at  the  higher  angles.  One  primary  cause  of 
this  is  the  crossflov;  lift  component  which  is  due  to  lateral  flow  separa¬ 
tion  and  the  formation  of  free  vortices  on  the  upper  surface.  This  non¬ 
linear  dependence, is  analyzed  in  References  1  and  5  s^d  summarized  by  Eaton 
in  Reference  6, 

The  expression  for  the  total  wing  lift  based  on  an  arbitrary 
reference  area  is 


+  C  +  C 
^BW-a 


(2) 


In  order  to  provide  high  angle  of  attack  capability,  it  is  necessary  to 
modify  both  the  linear  and  nonlinear  theories.  The  lift  on  the  wing  in 
the  presence  of  the  body,  as  presented  in  Reference  2,  is  a  linear  function 
of  angle  of  attack  and  can  be  expressed  as 


(3) 


Since  the  lift  force  does  not  vary  linearly  with  angle  of  attack  at  high 
angles.  Equation  (3)  is  modified  such  that  the  linear  lift  becomes  a 
function  of  sin  a  as  shown  below 


WB-a 


W 


(h) 
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This  ccjnponent  of  the  linear  lift  is  modified  further  to  satisfy  the  end 
condition  of  zero  lift  at  ,90  degrees  angle  of  attack.  The  resulting 
expression  for  the  linear  wing  lift  in  the  pressure  of  a  body  is 


C 


^WB-=a 


w 


sin  a  CQS  a 


(5) 


Ir  is  important  to  note  that  for  small  angles  of  attack  the  modified 
theory  should  he  very  close  to  the  method  of  Reference  1  since  sin 
a=a  and  cos  d^l,  Similarly,  the  additional  lift  on  the  body  due  to 
the  presance  of  the  wing  is 


BW-a 


*W 


sin  a  cos  a 


(6) 


The  parameters,  and  are  determined  from  Figure  5.  is 

obtained  from  Figure  6  by  multiplying  by  the  aspect-ratio 

ar"'*" 

if  AR  <  1.0.  When  the  aspect  ratio  is  greater  than  one,  the  lift-curve 
slope  is  obtained  from  the  following  ecjuation 


.  The  first  term  of  Equation  (?) 

is  an  empirical  modification  of  the  lift  curve  slope  for  a  lifting 
surface  with  aspect-ratio  greater  than  1, 


where 


is  '-'btained  from  Figure  ' 


k 


The  nonlinear  wing  lift  from  Reference  6  is 


=  C,  sin^a  /s„  \ 

hiY  ^ 

VREF/ 


(8) 


This  expression  is  modified  to  satisfy  the  aforementioned  end  condition 
with  the  result  being 


(9) 


where  is  obtained  from  Figure  To  It  should  be  noted  that  the  cross- 
flow  drag  coefficient  is  not  appreciably  affected  by  Mach  number, 
(References  6  and  7)»  hence,  is  presented  independent  of  Mach  ntombero 
The  total  tail  lift  is  computed  basically  the  same  as  the  wing  lift 
except  for  the  lift  due  to  deflection  of  the  control  surfaces »  The  tail 
lift  is  expressed  as 


C  =  C  +  C  +  C  +  C  +  C  (10)’ 

T  TB-a  BT-a  TB-5  BT-6  TV 


where  and  C  are  obtained  by  applying  Equations  (5)  and  (6) 

^TB-a  ^BT-a 

to  the  tail  surface.  The  other  three  components  are: 


'L  „  ^B  sin  6  f^T  \  cos{qi  +  6) 

'^ref; 


^T- 


sin  cos(a  +  6) 


(11) 

(12) 


Ct  =  Cj  sin^(a  +  6)  /s^  \  cos (a  +  6) 

^REF/ 


(13) 
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The  parameters,  eind  ,  are  obtained  from  Figure  5»  while  and 

are  obtained  as  specified  for  the  wing.  Notice  that  the  nonlinear 
L 

“t 

lift  is  based  on  the  local  angle  of  attack,  (a  +  6),  of  the  control 
surface , 

The  lift-loss  on  the  aft  surface  due  to  downwash  from  the  forward 
surface  is  obtained  from  the  method  presented  in  Reference  1  and  discussed 
in  Reference  ,  Since  the  method  for  computing  this  component  of  the 
total  lift  on  the  missile  is  both  complex  and  lengthy,  only  the  equations 
necessary  to  compute  the  lift-loss  are  presented.  The  reader  is  referred 
to  Reference  1  for  a  detailed  discussion  of  the  assumptions  and  technique 
used  in  deriving  the  method.  It  is  noted  that  the  nomenclature  used  to 
describe  the  lifting  surfaces  is  changed  from  wing  and  tail  to  forward  and 
aft  surfaces.  This  is  necessary  because  the  control  surface,  whether  it 
be  wing,  canard,  or  tail  type  of  control,  is  designated  the  tail  and 
because  the  aft  surface,  regardless  of  the  mode  of  control,  is  the  one 
which  is  affected  by  downwash. 

This  method  is  valid  for  the  entire  speed  range.  The  lift-loss  due 
to  downwash  is 


Kpg  sin  a  +  sin  6^, 

i(b  -  r)^ 

=  “f  “a 

- 

(iM 


This  equation  is  obtained  from  line-vortex  theory  assuming  only  one 
trailing  vortex  per  forward  panel  exists  (see  Figure  8),  The  lateral 
location,  f"  ®  and  the  vertical  location,  h.,  of  the  vortex  are  required 
for  use  in  Equation  (l4)  and  to  compute  the  interference  factor,  i. 

The  lateral  location  of  the  vortex  on  the  forward  surface  expressed  as 
a  fraction  of  the  exposed  semispan  of  this  surface  is 
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For  convenience,  the  right  hand  side  of  this  equation  is  defined  as  A', 
Isolating  f^,  results  in  the  following  expression 

fp  =  A'(b  -  r)p  +  Tp  (l6) 

where  f„  is  the  spanvise  location  of  the  vortex  at  the  forward 
£ 

surface.  Since  the  lateral  location  of  the  vortex  with  respect  to 
the  body  axis  is  unchanged,  the  subscript  may  be  dropped  and  Equation 
(l6)  may  be  written  as 


f  =  A'(b  -  r)p  +  rp  (IT) 

The  vertical  location  of  the  vortex,  h^,  is  measured  normal 
to  the  body  axis  at  the  center  of  pressure  of  the  aft  surface. 

The  expression  for  h^  is 

•“a  °  -'“r  -  W  S  [^A  *  <^CP>A  - 

Note  that  the  vertical  location  of  the  vortex  is  a  function  of  both 
angle  of  attack  and  the  deflection  angle  of  the  forward  surface. 

The  interference  f 8,ctor ,  i ,  is  given  by 


and 


f. 

X 


fr^ 


h.  = 
1 


\  t 
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Hence,  once  the  location  of  the  vortex  is  determined  and  the  position  of 
the  image  vortices,  f^  and  h^,  is  determined  from  Equation  (21),  the 
interference  factor  can  he  computed  from  Equations  (19)  and  (20),  The 
lift-loss  due  to  interference  can  then  he  calculated. 

The  total  lift  based  on  an  arbitrary  reference  area  is  obtained 
by  adding  the  components 


C 


L 


(22) 


where  ,  a  lift-loss,  will  be  a  negative  quantity,  The  other  three 

i 

components  include  both  linear  and :  nonlinei2,P  cpnt^ibuticms . 

DRAG  CHARACTERISTICS 

The  total  aerodynamic  drag  acting  on  a  missile  is  the  sum  of  the 
zero-lift  drag,  the  induced  drag  due  to  angle  of  attack  and/or  control 
surface  deflection,  and  the  base  pressure  drag.  It  is  well-known  that 
the  selection  of  the  proper  technique  for  computing  the  zero-lift  drag 
is  determined  by  the  operating  speed  of  the  missile.  Hence,  the 
methods  employed  to  compute  the  zero-lift  drag  of  a  missile  are  described 
for  three  speed  regimes.  These  speed  regimes  and  their  associated 
limits  are  defined  as  follows ; 


1,  Subsonic  —  M  <  0,8 

2,  Transonic  —  0,8  -  M  -  1.2 

3,  Supersonic  —  M  >  1.2 

The  description  of  C^^  calculations  for  these  speed  regimes  is  followed 

o 

by  a  description  of  the  method  used,  to  compute  the  induced  drag  due  to 
angle  of  attack  and/or  tail  deflection.  The  last  section  presents 
the  total  drag. 
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lo  Subsonic  Region 


The  zero-lift  drag  for  subsonic  speeds  can  be  expressed  as 


(23) 


where  each  of  the  components  is  composed  of  skin- friction  drag  and  pressure- 
drago  The  pressure  drag  at  subsonic  speeds  is  usually  small  comparedto  the 
drag  due  to  skin  friction,,  Since  the  flow  around  high  speed  missiles  of  the 
type  being  considered  here  is  primarily  turbulent,  the  methods  employed  are 
developed  assuming  the  existence  of  fully  turbulent  boundary  layers o 

The  zero-lift  drag  of  the  body  based  on  an  arbitrary  reference  area  is 
obtained  from  Section  ko2o3<.l  of  Reference  9  and  can  be  expressed  as 


C  =  1„02  C  |1  + 


+  7 


REF 


(2k) 


where  C^,  the  skin- friction  coefficient,  is  determined  using  Figure  9= 

The  wing  zero-lift  drag  as  presented  in  Section  4olo5ol  of  Reference  9 
is 


8o0  C. 


w 


w 


1  +  2(t/c)  +  100(t/c) 


w 

""ref 


(25) 


The  factor  of  8.0  is  included  to  accoimt  for  the  total  wetted  area  of 

the  surfaces  and  for  the  existence  of  four  wings.  The  tail  zero-lift 

C. 


drag. 


D 


can  be  obtained  by  using  the  tail  thickness-to-chord  ratio  and 


the  total  tail  area,  S^  ,  in  Equation  (25). 

These  three  components,  based  on  the  same  arbitrary  reference  area, 
are  added  as  indicated  by  Equation  (23)  to  give  the  total  subsonic  zero- 
lift  drag. 
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2,  Transonic  Region 

The  "body  drag  in  the  transonic  speed  regime  is  ccsnposed  of  compressible 
skin-friction  drag,  subsonic  pressure  drag,  and  transonic  wave  drag  and  is 
obtained  using  the  techniques  presented  in  Section  of  Reference  9o 

The  compressible  skin- friction  drag  is  obtained  using  the  following 
equation: 


=  1.02  C-  S 
f  s 

REF 


(26) 


where  is  a  function  of  both  Reynolds  nxanber  and  Mach  number  and  can 

c 

be  determined  from  Figures  9  an<i  10.  The  subsonic  pressure  drag  as  ex¬ 
tracted  from  Equation  (24)  is 


=  1.02  C 

B 


1,.5 


s 

s 

®REF 


(27) 


The  transonic  pressure  drag  is  computed  using  Equation  (27)  up  to  a 
Mach  number  of  1.0  and  then  it  is  decreased  linearly  from  its  value  at  1.0 
to  zero  at  M=lo2o 

The  transonic  wave  drag  of  the  body  is  obtained  from  Figure  11  which 
presents  the  wave  drag  as  a  function  of  the  nose  fineness  ratio  and  Mach 
number.  This  figure  was  constructed  from  experimental  data  presented  in 
Reference  10  and  by  using  the  curves  of  wave  drag  for  bodies  of  revolution 
as  shown  in  Reference  11  and  reproduced  here  in  Figure  12. 

The  total  transonic  body  zero-lift  drag  is  obtained  form  the  following 
expression: 


+  C. 


+  c.. 


B 


V  S 
B  REF 


ib 


(28) 


Experimental  results  show  little  increase  in  the  viscous  drag  of  the 
aerodynamic  surfaces  from  the  subsonic  to  the  transonic  regime  and  there¬ 
fore,  the  skin- friction  drag  for  the  subsonic  region  is  also  used  in  the 
transonic  regime o  It  may  be  expressed  as  follows; 


(29) 


To  this  wing  transonic  skin- friction  drag  is  added  a  drag  increment ,  AC^  , 

°W 

which  is  the  transonic  wave  drag  of  the  wing  surfaces„  Figure. 13  expresses 

this  component  as  a  fxmction  of  thickness-to-chord  ratio,  —  ,  aspect-ratio, 

c 

and  Mach  number  for  rectangular  surfaces „  For  surfaces  having  swept 
leading  edges ,  is  obtained  as  for  rectangular  surfaces  and  then  ad- 

°W 

justed  to  account  for  the  sweep  angle  using  the  following  equation; 


AC 


D 


cos 


(30) 


The  Mach  number  jus ed  in -FJ.gure  13  to  obtain  A  for  swept  lifting 
surfaces  is 


M  =  M 


cos 


^c/uj 


1/2 


(31) 


This  component  of  the  drag  due  Ji-O’  the  ass'rodynamic  surfaces  must  be  com¬ 
puted  for  both  the  wings  and  tails,,  The  tail  contribution  is  determined 
in  the  same  manner  as  the  wing  contribution  above  using  tail  parameters,, 
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The  total  transonic  zero-lift  drag  is  the  sum  of  these,- components  as 
shown  below: 


3.  Supersonic  Region 

A  simple  empirical  method  for  computing  the  zero-lift  drag  for  the 

supersonic  speed  regime  has  been  developed  by  assuming  a  parabolic 

variation  of  with  Mach  number  between  1.2  and  3.0.  The  resulting 
o 

equation  which  is  used  to  compute  the  zero-lift  drag  of  a  missile  for 
Mach  numbers  greater  than  1,2  is 


p 


1  ~  (3/ 


(33) 


where  '  and  "  are  the  values  of  the  total  zero-lift  drag  at 
^o  o 

Mach  numbers  1.2  and  3.0  respectively.  It  should  be  noted  that 

although  for  Mach  nmbers  greater  than  3«0  can  be  determined  from 
o 

Equation  (33),  existing  hypersonic  flow  theories  would  probably  provide 

a  more  accurate  estimate  of  the  zero  lift  drag. 

In  order  to  utilize  Equation  (33)  for  determining  the  variation 

of  with  Mach  number,  '  and  "  must  be  specified,  '  is 

o  o  o  o 

determined  by  using  the  techniques  described  in  the  previous  section 

for  the  transonic  flow  regime.  Since  the  magnitude  of  the  supersonic 

wave  drag  is  heavily  dependent  an.. the  nose  shape  of  the  missile,  " 

is  determined  using  one  of  two  methods;  the  selection  of  the  proper® 

method  depends  on  the  missile  forehody  shape. 
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In  the  first  method,  which  is  for  blunted  ogives,  pointed  ogives, 
and  blunted  cones,  the  supersonic  zero-lift  drag  is  a  fiinction  of  Mach 

II 

number  and  nose  fineness  ratiOo  Since  is  the  zero-lift  drag  at  M=3o0, 

it  remains  to  define  its  variation  with  °nose  fineness  ratiOo  Hoerner, 
Reference  12,  .indicates  the  zero-lift  drag  for  body-fin  configurations  with 
slender  (high  fineness  ratio)  nose  shapes  generally  peaks  at  a  Mach  n-umber 
of  loO  to  lo2  and  then  decreases  to  approximately  its  subsonic  value  plus 
the  transonic  wave  drag  of  the  lifting  surfaces  at  M=3„0o  Similar  con¬ 
figurations  with  blunted  nose  shapes  of  low  fineness  ratio  reach  a  peak  at 
about  the  same  Mach  number,  but  decrease  very  little  as  the  Mach  nimiber  is 
increasedo  Using  these  trends  as  an  indication  of  the  effect  of  fineness 


ratio  on  the  variation  of  zero-lift  drag  with  Mach  number  for  the  supersonic 

II 

speed  regime,. for  the  aforementioned  nose  shapes  is  specified  as 
follows ;  ° 


a)  For  (l/d)j^  <  0„5, 

b)  For  (l/d)jj  >  8o0, 


where 

o 


is  determined  utilizing  subsonic  flow  theory  and  A 


ri 


and  ACj^  are  the  transonic  wing  and  tail  wave  drago  The  variation  of  the 
°T 

forebody  wave  drag  as  a  function  of  Mach  nmber  is  presented  in  Figure  12  and 
was  used  to  construct  Figure  l4  which  is-utilized  with  the  following  equations 


(34) 


tr 

to  compute  for  nose  fineness  ratios  between  0o5  and  8o0o 


In  the  second  method,  which  is  for  pointed  conical  noses,  C  is 

o 

determined  from  the  following  equation; 


c  =  c  / 

%  ^o/m=0o8 


+  /  AC. 


W  +  AC 


°W  ^REF 


°T  ^REF/  Tr 


where  the  first  two  terms  are  obtained  in  the  same  manner  described  above. 
The  forebody  wave  drag,  which  for  this  type  of  nose  is  a  function  of  Mach 
nxmiber  and  the  cone  seraivertex  angle,  is  obtained  from  Figure  15. 

If 

Once  Cjj  is  determined.  Equation  (30)  can  be.  used  to  compute  the 
o 

zero-lift  drag  at  any  Mach  number  between  1.2  and  3.0.  The  complete  zero- 
lift  drag  curve  of  a  missile  can  now  be  determined. 

Uo  Induced  Drag 

The  induced  drag  due  to  angle  of  attack  and/or  tail  surface  de¬ 
flection  is  composed  of  four  drag  increments  as  shown  below: 


The  drag  increment  for  the  body  is  obtained  using  a  method 'presented  in 
Reference  3,  while  the  induced  drag  due  to  the  wing  and  tail  is  obtained 
from  the  drag  of  an  equivalent  flat  plate  normal  to  the  flow.  The  in¬ 
duced  drag  on  a  body  of  revolution  at  angle  of  attack  can  be  expressed 
as  follows ; 


The  wing  induced  drag  "based  on  an  arbitrary  reference  area  is  obtained 
from 


(38) 


where  is  the  drag  of  a  flat  plate  normal  to' the  flow  field,  and  is 

^FP 

determined  from  Figure  l6o  The  drag  curve  shown  in  Figure  l6  has  been 
constructed  using  the  three-dimensional  subsonic  drag  coefficient  for  a 
flat  plate  normal  to  the  flow  (Reference  13)  and  the  variation  of  the  drag 
coefficient  with  Mach  number  for  the  two-dimensional  flat  plate  (Refer¬ 
ence  12)  „  The  equivalent  flat  plate  area  of  the  wing  is  taken  as  the  pro¬ 
jection  of  the  wing  area  on  the  normal  plane,  sin  Similarly,  the 
drag  increment  of  the  tail  surface  at  angle  of  attack  is 


AC 


T-a  ^ 


sin  a 


REF 


(39) 


The  drag  increment  due  to  deflection  of  the  tail  surface  cannot  be 
obtained  directly  since  the  drag  increment  is  not  a  linear  function  of 
local  angle  of  attack,,  The  total  drag  increment  for  the  tail  surface  can 
be  expressed  as 


^T-a  "  ■^T-6  ^FP 


'sin(a  f  5) 

'■"ref 


'(Uo) 


where  (a  +  5)  is  the  local  tail,  angle  bf  attacko  It  now  becomes  a 
simple  matter  to  obtain  the  drag  increment  due  to  the  deflection  of  the 
control  surface  by  taking  the  difference  between  Equations  (39)  and  (U0)o 
Thus 


AC„ 


T-6 


T-a 


(1+1) 
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6.  Total  Drag  ■ 

The  total  drag  can  be  expressed  as  follows  J 

S  ^  %.  (42) 

where  is  determined  using  Equation  (36)  and  is  obtained  from 
i  ^o 

Equation  (23) »  (32),  or  (33)  depending  on  the  speed  of  the  missile. 

Base  pressure  drag  is  not  included  in  Equation  (42), 


l6 


PITCHING  MOMENT  CHARACTERISTICS 

The  total  pitching  moment  acting  on  the  missile  is  the  sum  of  the 
moments  due  to  the  lift  and  drag  forces  acting  on  the  body,  wings,  and 
tails o  Most  methods  for  computing  the  pitching  moment  (References  1,  6, 
and  9)  consider  only  the  moment  due  to  lifto.  This  is  valid  only  for 
small  angles  of  attacko  If  large  angles,  of  attack  are  to  be  considered, 
the  moment  must  include  the  drag  contribution,,  In  general,  the  body 
longitudinal  pitching  moment  is  determined  directly;  while  the  other  com¬ 
ponents  are  determined  only  after  the  centers  of  pressure  of  the  wing 
and  tail  surfaces  are  specified„ 


The  body-alone  pitching  moment  about  its  center  of  gravity  is 
obtained  from  the  method  of  Allen,  References  2  and  3o  The  expression  is 


where  (k^  -  k^) ,  R  determined  from  Figures  3  and  This 

moment  coefficient,  C  ,  is  based  on  an  arbitrary  reference  length  and 

• 

are  So 

As  noted  above,  the  center  of  pressure  locations  for  the  wing  and 
tail  surfaces  must  be  specified  before  their  pitching  moments  can  be  de- 
terrainedo  It  must  be  remembered  that  the  linear  lift  of  the  aerodynamic 
surfaces  is  composed  of  two  components;  'the  lift  on  the  surface  in  the 
presence  of  the  missile  body  and  the  additional  lift  on  the  body  due  to 
the  presence  of  a  lifting  surface „  This  means  that  in  order  to  deter¬ 
mine  the  linear  pitching  moment  caused  by  the  lifting  surfaces,  it  is 
necessary  to  specify  the  center  of  pressure  location  for  each  of  these 
linear  components o 

The  center  of  pressure  of  the  lift  on  the  wing  in  the  presence  of 
the  body  as  measured  from  the  junction  of  the  wing  leading  edge  and  the 
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body  is  obtained  using  Figure  17 o  The  reference  point  is  transferred  to 
the  nose  by  using  the  following  equation 


where 


is  obtained  from  Figure  17a  for  subsonic  speeds  and  from 


Figure  17h  for  supersonic  speeds. 

The  center  of  pressure  of  the  additional  lift  on  the  body  in  the 
presence  of  the  wing  is  obtained  using  Figure  l8  if  the  flow  is  subsonic, 
and  either  Figure  I9  or  20  if  the  flow  is  supersonic.  For  the  case  of 
supersonic  flow,  Figure  19  is  used  if 


3AR  (1 


\ 


and  Figure  20  is  used  if  the  above  quantity  is  greater  than  U.O.  The  cen¬ 
ter  of  pressure,  as  obtained  from  the  aforementioned  figures,  is 

referred  to  the  nose  by  using  the  following  equation. 


.  .The  location  of  both  centers  of  pressure  for  the  tails,  and 

(^Cp)g^,  can  also  be  obtained  from  the  above  procedure. 

The  centers  of  pressure  for '-a  given  TtTting  surface  are  combined  to 
obtain  a  single  average  center  of  pressure  location  for  each  set  of  aero¬ 
dynamic  surfaces.  For  example,  the  average  center  of  pressure  of  the 
wings  is  obtained  by  computing  the  total  pitching  moment  due  to  the  wings 
and  dividing  by  the  wing  normal  force,  -‘The  pitching  moment  about  the 


tils, 
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nose  of  the  body  due  to  the  wing  is 


\  ^  H/v  h) 


C  cos  a 

Sw 


AC^  sin  a 
^W-ct 


(W) 


where  all  of  the  above  terms  have  been  previously  definedo  The  average 
wing  center  of  pressure  as  measured  from  the  nose  of  the  body  can  now  be 
defined  as 


(ii7) 


(  \  ^mW  ^REF 

I  Cp4  -  J  cos  a  . 


+  AC  sin  a 
^W-a 


Similarly,  the  tail  pitching  moment  and  center  of  pressure  including  the 
effect  of  surface  deflection  can  be  expressed  as 


^TB-6 


I  f  C,.  +  Ct  +  C  \  cos  a 

[v  ^B-a  ^TV  .  ^ij 

(=‘cp)B/bE: 


:  cos  a  (1|8) 

L  ^BT-a 


where 


“^mT  ^REF 


and 


\  T-a  1/  .T~-6  \  „T-a  T-6/ 


c  =  c  +  c  +  c 

T-a  TB-a  BT-a  TV 


C  =  C  +  C 

T-6  ^TB-6  BT-6 


(i*9) 
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The  total  longitudinal  pitching  moment  of  the  missile  about  the 
missile's  center  of  gravity  may  be  expressed  as 


(50) 


COMPUTER  PROGRAM  DESCRIPTION 

The  method  presented  herein  Tor.  obtaining  the  static  aerodynamic 
characteristics  of  a  missile  has  been  programmed  for  use  on  the  IBM  7090 
computer  and  other  compatible  digital  computers 0  The  program,  Table  1, 
is  written  in  Fortran  II  and  requires  only  the  geometric  characteristics 
of  the  missile,  and  its  flight  conditions  as  inputs.  The  output  consists 
of  the  static  longitudinal  aerodynamic  characteristics  in  coefficient 
form,  the  center-of-pressure  location  for  the  body,  wings,  and  tails, 
the  lift-curve  slope  for  each  set  of  lifting  surfaces,  and  the  compo¬ 
nents  of  the  lift  and  normal  force  coefficients.  The  force  components 
referred  to  are  the  body,  wing,  and  tail  lift  and  normal  force  coeffi¬ 
cients  and  the  coefficient  representing  the  lift-loss  due  to  downwash. 

The  aerodynamic  characteristics  are  output  in  both  the  stability  and 
body  axis  systems  (Figure  21).  Provision  has  been  made  for  a  third  lift¬ 
ing  surface  to  account  for  the  possibility  of  using  strakes  in  combina¬ 
tion  with  two  other  sets  of  lifting  surfaces  (Figure  l). 

The  program  itself  consists  of  a  main  program  and  three  subroutines- — 
GEOSUB,  CLASUB,  and  CATSUB.  The  first  subroutine  performs  some  initial 
geometric  computations,  determines  the  nose  wave  drag  constant  (Figure  lU), 
and  obtains  the  Reynolds  number  per  foot  based  on  the  altitude  input  to 
the  program.  Subroutine  CLASUB  determines  the  lift-curve  slope  of  the 
lifting  surfaces  from  curve  fits  employed  to  represent  the  curves  pre¬ 
sented  in  Figure  6.  The  last  subroutine,  CATSUB,  obtains  the  body-wing 
and  body-tail  interference  factors,  computes  the  center  of  pressure  loca¬ 
tion  as  a  function  of  the  root  chord  of  the  lifting  surfaces,  and  deter¬ 
mines  the  crossflow  drag  coefficient  for  the  lifting  surfaces. 

The  program  computes  the  static  force  and  moment  coefficients  for 
typical  missile  configurations  at  specified  angles  of  attack,  control  sur¬ 
face  deflection  angles  and  Mach  numbers.  The  angle  of  attack  range  is 
-180°  to  +180®,  and  any  control  surface  deflection  within  this  range  may 
be  used.  The  Mach  number  is  limited  to  3.0  only  because  the  drag  pre¬ 
diction  methods  are  valid  up  to  this  particular  Mach  number.  The  pro¬ 
gram  can  be  used  for  configurations  at  M  >  3.0;  however,  the  drag,  pre¬ 
dictions  above  this  limit  should  be  used  with  caution.  The  computer 
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program  can  "be  used  to  obtain  build-up  information;  that  is,  the  aero¬ 
dynamic  characteristics  of  the  missile  body  alone,  the  body-wing  configu¬ 
ration,  and  the  body-tail  configuration.  This  information  may  be  ob¬ 
tained  by  simply  setting  the  appropriate  parameters  to  zero. 

The  inputs  to  the  computer  program  with  their  FORTRAN  symbols  are 
presented  in  Table  2,  The  format  for  preparing  the  input  cards  is  pre¬ 
sented  in  Table  3o  It  should  be  noted  that  if  a  configuration  does  not 
have  a  control  surface,  e.g,  a  body  alone  configuration,  the  number  of 
control  surface  deflection  angles  should  be  set  at  one  and  the  deflec¬ 
tion  angle,  itself,  would  be  0,0  degrees.  Cards  11  and  12  are  used 
only  when  the  number  of  angles  of  attack  require  their  use.  The  output 
variables  are  defined  in  Table  4,  There  is  no  limit  to  the  number  of 
data  decks  which  may  be  stacked  together  and  run  at  the  same  time. 


COMPARISON  OF  THEORY  WITH  EXPERIMENTAL  DATA 

Mumerous  comparisons  of  theory  with  experiment  have  been  made  in 
order  bo  establish  and  verify  the  accuracy  of  the  method.  Figure  22 
presents  the  configurations  used  for  comparison.  Configuration  1  is  a 
strake-tail  configuration.  Figure  23  presents  the  comparison  between 
the  experimental  data  and  the  theoretical  results  obtained  from  the 
method  described  herein.  The  theoretical  lift  and  drag  are  within 
15  percent  of  the  experimental  data  for  M  =  0,T»  The  deviation  for 
M  =  1.1  does  increase  to  about  25  percent.  However,  it  should  be  noted 
that  the  accuracy  of  the  method  to  30  degrees  is  good.  The  method  also 
predicts  the  center  of  pressure  location  very  satisfactorily. 

Configuration  2  is  a  wing-controlled  vehicle.  Comparisons  are 
presented  in  Figure  24  for  Mach  numbers  of  1.12  and  2.l6  and  for  wing 
deflection  angles  of  0  and  -10  degrees.  Although  high  angle  of  attack 
data  is  not  available  for  this  configuration,  it  appears  that  the  com¬ 
paratively  high  theoretical  lift  shown  by  Configuration  1  may  be  more 
pronounced  for  Configuration  2,  This  trend  is  seen  by  the  decrease  in 
the  slope  of  the  experimental  lift  data  at  25  to  30  degrees  angle  of 
attack,  Figure  24.  The  comparison  with  the  wing  deflected  down  10  de¬ 
grees  appeared  to  be  good. 
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Data  for  Configuration  3  is  presented  in  Figure  25 ^  Configuration 
4  was  tested  dy  the  Cornell  Aeronautical  Laboratory  to  180  degrees  angle 
of  attack.  Reference  l4.  Comparison  of  theory  with  experimental  data 
is  presented  in  Figure  260 


CONCLUSIONS 

A  method  for  predicting  the  static,  longitudinal  aerodynamic 
characteristics  of  low  aspect-ratio  missiles  operating  at  angles  of 
attack  to  I80  degrees  has  been  developed^  The  method  is  valid  for 
a  wide  speed  range  and  considers  control  surface  deflections c  A  com¬ 
puter  program,  written  to  facilitate  use  of  the  method,  has  been 
described,,  Results  obtained  using  the  method  have  been  compared  with 
wind  tunnel  data  and  acceptable  agreement  has  been  demonstrated^ 
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Figure  1  —  Typical  Missile  Configurations 
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Figure  5c  —  (^AR)(1  +  A)(l/m/3  +  1)  >  4;  No  Afterbody 
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Figure  6  —  Lift  Curve  Slope  for  Wings  and  Tails  (from  Reference  8) 
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Figure  6b  —  Unswept  Mid  Chord 
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Figure  13  ;—  Transonic  Zero-Lift  Wing  Wave  Drag  for  Unswept  Wings  (from  Reference  9) 
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Figure  18b  —  No  Trailing  Edge  Sweep 


li8 
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Figure  19a  —  No  Mid  Chord  Sweep 
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Figure  I9b  —  No  Trailing  Edge  Sweep 
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i  Figure  21  —  Missile' Axis  Systems 
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CONFIGURATION  4 

Figure  22  —  Configurations  Used  to  Compare  Theory  with  Experiment 
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Figure  2,3  —  Comparison  of  Experimental  Data  with  Theoretical  Results  for  Config 
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Figure  23  (Continued) 
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Figure  24  (Continued) 
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Figure  24  (Continued) 
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Figure  25  (Continued) 
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Figure  26  —  Comparison  of  Experimental  Data  with  Theoretical  Results  for  Configuration  4 


Table  1 


COMPUTER  PROGRAM  LISTING 


Dlf^ENSleN  XVXMrl63,XDTri6]^XALC483 
XVxM#XDTiXAL 

ce^M'^N  CNiCA^CNB,CNWiCNTiC^-.'W2/f:LW/CLT/Cl.W24CLl#CLWB/CLVlSW,CLlT/ 

1  CUlWiCNTD,CABiXCP2iyCG2/XCPBiXCPT^XCPW/XCPW2/XCG/XCPTV/ 

a  xlamta 

Cer^MOM  LLKKi  iSwPW,  I  SWPT/  I  S'*'PW2i  IAFBWi  I  AFbTj  I  AF8W2>  IL^LLLL#  IJ/ J/ 

1  I  I^NBBDY#  IZZYi  ICSC^  iNOSE,N:M#Nl^LK/  IDT#  I  AL/ISWPl  /  I  AFB 

CeMM0|Sj  xLaYI^'#Xl  AMT^  XLA^W2>yMAC'*''/XMAcT  /  XMaCW2^CLAMW|CUAY!T^  Ci.AY'Wpi 

1  EW,BT*Bu2^CR0eTW>CRBOTTJ>CP9f’w2/SW,ST/SW?^yWrNG/XTAlt.,  XWING2 

2  /XL^DiDliXLN'/ARFA/XPFP^SsUBS^XLBBiZFiART^  ARW/ARW2/XLSj9SE' 
C8-^0N  CBLAYI/BCPLAM/CReeT/Pl/CLALUxlA^l^BARiRATie^XKTB/XKBT/ 

1  XKwB/XKgW^  XBCRB'-''i  yKWBW^  XkWBT/  XKWBW2<  XKBWW*  XKBWT#  XKBWW2/ 

2  XCpBW^ XcPBT,XCPBW2^ vCPWB,XCPtS#^CPWB2>XKWBI/8CC^ ODCW/flDCTj 

3  ePcW2^CLALW^  CLAL.T/CLALW2^REFT/BFTA/ AU/  TONSTjiHBNiST^  PT/  XKTBI  / 

4  VXM/VXMRl,D2LTA/XKc'?WB^XBCPWQiXHAC>XLAMW4/XLAM24/XLAM2/ 

5  XLA^A/  TBVC/  revCT/  TR\/CW2/  TBVCW/  EXS#  STT0T#SWTBT<  S'v'PTBT/  RE/ 

6  CPR>CDe'rt/CD0T^CP8B/rDftV«2^CL8W^CLWB2>CLBW2/CLTRiCLBT^CLTD^ 

7  CLTDeiCLBDT^CLVlST>rLVlW24CMB/CD0WBT 
reMMeN  DEoBsw^nCDesT^DCDosP/ncoeWiBcDeT^cCDewp 
ceYMBN  cdaUZ 

lU  '  0 

3C4C  FSRf^AT  [2I5>  7F1c«6] 

3C5C  PePR*TC6l53 
3C2C  FRRYAt C7P10«3] 

3C1C  FeRY:AT  CICaA] 

311C  FeP^AT  [IHI/ 10A63  ,,  ^ 

312C  F0RYAT C6X,2HHTy  9X/2MD  / 8X^ 2HXL, 6X/ 6HXLN9SE/ 5X, 3WXCG/ 6X/ AHaREA / 6Xi 
14RXREF] 

3C21  PSRYA  j  [//,  5X/  5MT9VCWi5X/  6l-iT0VCW2^5Xi5HTPvCT3 
314C  F9RRAT  C  3X^  l5,  5X^  15,  kx*  7F15'6#//'] 

315C  RgRf^AT  C4X,  I5,5X/ I5,5X,  IB/SXi  I5^6X/ I5/5X/ 15] 

324C  FeR^AT  C  6X/5mISWPW/  5X,cHIaFbW,  lOX#  5HXLAMW,  IQX,  5MCLAMW, 

1  1CX,5H  Ew  ^  9v/6HCRB3m/ KX,58  SW  /  ICX,  BHX^ACW  ,  1  OXi  5MXW  I  NG^ /] 
3244  fSRRAT  C  6Xi6FlSWPW2,  4X,fMIAPBW2^  9X/  f>HXLAMW2jt  9X,  feHC.tAMW2, 

1  1CX,5M  Bw2  ^  9X>6HCR80'n2/ 10X,5H  5W2  ^  9X,  6HXMACW2, 9X#  6HXWINg2/ /] 
32Ag  FBRf^AT  [  6X,5HISwPT/  BX/BHIAPBT/  5hXLAMT^  IQX/  5HCLA''‘T/ 

1  IOX45P  Bt  / 9Xi 6HCR03TT/ IFX, 5H  ST  1  lOX,  BPXMACT^ lOX# 5HXT A  I L/ /] 

3241  reRMAj  C//,6Xi5PlCSC  ,5X/5HTN0SE,5X^5PNDELTj5Xi5MNMACM,5X/5HNALPH, 

IBXiBPNOeDY]  ■  „ 

3333  read  3010,  T  I  Tl  1 ,  T I  TL2',T  I  Tl3,  T  I  TLA  /  T  I  TL.5,  T  I  TU6/ T  T  TL7i  T  I  TLSi  T  I 
ITITLC' 

read  305c,  ICSC/ IN0SE, IDT, IM, IAl,NBODY 

read  3C4C,  ISWPW,  I  AFPW,XLA"W,CLAMW,Bl^'/CR00TW,SW,XMACW,XVnNG 
read  304  0,  ISWpWp,  I  AFB''’2,XLAT'W2,CUAmW2/BW2,CR00W?,SW2/XMACW2,XWIN 

1G2 

RFAD  304C,  ISWPT,  I  AFBT/XLA‘'’T,CLAMT4BT/CR6eTT>ST,XMACT,XTAlL 
read  3020,  HT,C/XL,XLlSi0SE,yCG/AREA,xPEP 
read  3020,  T0VcW,TevCw2/TevCT 
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PRI^'-'T  3llo/ 
ITITLO 
print  32'+1 
print  3l5o/ 
print  3240 
print  3140/ 
print  3344 
print  3140/ 
1XWING2 
print  3248 
print  3140/ 
print  312C 
PRINT  3C2C/ 
print  3021 
print  3C2C/ 
IL  *  1  + 


T ITLI/ T ITl2/ TiTL3/TiTU4/TITLe/TITL6/TrTL7/ TITUS/ TITL9/ 

ICSC/ IN0SE. IDT/ IM/ I AL/Nb8DY 

ISUPW/  I  AFBW/XL  A'^W/CLAHW,BW/CRe9TW/SW/XMACW/XWlNG 
ISwPW2,  IAFPW2,yLAMW2^CLa'^W2/BW2/CRB0W2/SW2/XNACW2/ 

ISKPT/IAFBT/XLAYT/CLAMT/BT/CRSe'TTiST/XMACT/XTAIL 

HT/D,XL/XLN‘eSE/XCG/ AREA/XREF 

TevCW/TevCN'2/TBVCT 

IL 


llkn=o 
llll=g 
XCGE*XcG 
IZZY  =  C 
caul  GFBSuB 

read  40C0,  [XDT CNl ^Nai, IDT3 
read  40C0#  CXVXN  [K] IMI 
read  4OCO,  CXALCNAI /NA*1/ IaL] 
P£bREFT-»VxN 


VXN’XVXN  Cl] 

00  60c2  Ij=l/  IN 
OEU'C AiexDt  Cl] 

08  6CC1  1I*1/I0T 
AlPHAbXAU Cl3 
4CCC  reRYAT C16f5. 1] 

5CCC  FaRY-AjriHl/4HVxYa/F5.2/2X^  aHDELTAbiF6*2//3 

5CC1  FgR^Aj  C2X^2HAL,3y/5HCLT0T^  pXj,5HcDT8t/2X/4HCLWT/3X/4HCLTT/3X/3HCUB/ 
14X/3WCLI/4X/ 3McNW/ 4X/ 3HCNT  »  4X/ AhCNTd/ 3X/ 3HCNB/ 5X#  2HCN/ 5X/ 2HCA/ 4X/ 
H4PXCPW/3X,4HXCPT,3X/4hXCPB; 3V/AhXCP2# AX/2HCM//] 


82CC:  print  ScCCi  VXN/DELTAl 
print  5CC1 

0ELTA=DeLtA1/57'29578+.C00O00001 
DP  6000  G«1/IAL; 

AL =AUpWA/57* 29578+. OCOGCOOl 

1  VX'^RIsVX't 
IZZY=IZZY+1 

IF  CiZZY  -  4]  6666/6666/1111 


6666  VXN*.6 
nil  CALL  CUASuB 

JFCLLLL-1]  900/942/980 
9CC  IFClZZY-4]  6009/6009/925. 

925  CALL  CATSUB 

6CC9  xLAN14*ATaN  CC.5^»  CBl  »C]  «  1 . /CSL AM^ , 25  +  XU AM1-»CR99T“  . 25-®CRB0T3  / 

1 C.5*Cb1-D]]] 
r£sREfT«VxM«XHaC 
IF  CRE-1.E06]  6OI 0/ 6020/ 6020 
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6C1C  AA»‘C835 
yNNs*.2ll 
G9  T8  6070 

6C2C  IF  CRE-I'EOTI  6030/6040^6040 
6030  AA»*05? 

yNN»-tl77 
oe  T0  6070 

6040  IF  CRE-1'EC8]  6050/ 6060/ 60a0 
6C5C  AA*'C333 

XNN*".1488 
GQ  Te  6070 
6060  AA»»0221 
y\N""tlE7 

6070  CF  =  ^A*RE:**yNN 

CCP''*  .«CF*  n  .■*-2**TpVC+lC0t*T8VC#«4.] 

rxs»  i:d*cO /  [8**ceLAM] 

IF  [IS'>^'P1»13  6c80/6080/6090 
6080  rXS*2,*E^S 

6090  EXS*  tCOEOT-*D/2.+FXs]  *2. 

IF  ClZZY-4]  6091/6051/6002 

6091  IF  ClZZY-3]  6093,6094/6055 

6092  IF  CAl]  2^^01/2^02/  2402 

2401  f'CC*-eOC 

2402  I.LXK  =  LLKK+1 

IF  CLL^X-El  24c3, 24C4, 2420 

2403  IF  CSw]  2410/2410/2420 

2410  LLKKclLKK+1 

2404  IF  CSw2]  2411/2411/2420 

2411  lLKK=LLXK+1 

2420  IF  CLLXK"2]  530/943/950 
930  vk-vB'WeXKWb 
X  K  B  W  E  y  K  B  w 

VCP^BeVy,  lNG  +  XBcRwB*fCR8PT 
vxPBWcXw  In,G  +  XBcRP.W»CR8RT 
PCC'**  =  9CC 

560  CL-‘‘  *S  I  f-j  C  alI  «  CXkWoW  +  XKB’^''W3  -ifrLALW^SW^cOS  Cal3  /  AREA 
CaUI  *XKWbW*CLALW*SW»C0S  [ALl  /  AREA 
f  L  B  ”  C  L  "  C  L  w  B 

CLVisw=  csiN  [AL] *s I N  caL3 »sw»ces  cal^ /aReaJ  «eDcw 

cl’**'  ®  Cl  •'^+cl'^  I  sw 

6093  CC9''"  =  C0e*  CSW  +  ExS] /area 

XLAf^WAsXLA^H 

TevcweTevc 

SKT9T=SW+EXS 

IZZY»IZZY+1 

IFlSWEl  S42/  94F/511. 

51 1  C6LAY  =  CPS[CLAMx2l/SIN  CCLAHl-'2] 
pCCLAM=BFTA*CSlAM. 
cR38T  =  r.Reew2 
9  1 5  B  W  2 
1AFB«IAFBw2 
CLALIeClAlWE 
XLA’^1  eXLAMW2 
Tq vC=t3VCw2 
XP  AC«x'^ACw2 
ISwPIeTSWpWB 
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O’)  rv 


eAR*BETA»ARW2 
fJATlfl»CRa0T/ CBETA«D] 

IFtIZZV-A]  6009^-6009/925 

943  yK'^BW2»XKwB  .  - 

XKB'aW2  =  XKBW 

ycP'A'B2  =  XWiNG2+xBcRWB«CR89T 
XCPB'*<2»xWiNG2+'XBCRDW#CR80T 
F'DCWB.acC 

944  cL’a2*SIN  CaLJ  •  CxKWBW2'>'XK6WW?]  «CU ALW2#SW2»c0S  CALD /AREA  ■ 
CLaE2»SIN[AU] *XKWBW2#cUALW2«SW2*C0ScAt3/AREA 
CLB'^2»CLW2-CLWe2 

CLVrw2»  CSINCAL]^^SI^CAL]»Sw?«C8S[AL3/AREA3^^eDCW2 
CL’'J2'=CLw2+CLVIw2 
6094  rC9w2rCDP*  [SW2+EXS]/AREA 
XLA^24=XlA^14 
VLA^2*XLAm14 
Sw2T«T=SW2+EXS 
S4E  IZZY  «  IzZY  +  1 
LLK<»LUKK+2 

IF  [ST]  90C/9gO,94O 
940  CSLA^^^CeS  [CLAYt] /SIK  [CLAMT] 

ART*  CbT-D]  ih»2/sT 

PC9LAf^«BETA^»C6EAM 

CROSTsCRPeTT 

di=bt 

eAR*BETA*ART 
CLALIsCLAlT 
I AFB= I AFBt 
XYACsx'^ACt 
TOVC=tPVCt 

ISa'FI*  iSWpT 
XLAf'UXLAMT 

RAT IP  =  CRP9T/ [BeTA*D] 

IF  [IzZY"A]  60C5# 6009/ 925 
950  XK'*vBT  =  Xk'a'b 
X  K  c  T  *  X  K  B  W 

xcPBT=XtAiL+XBcRBW*CR6PT 

"CCT=o9C 

551  CLT'  [  [XKWbT+XKeWt]  ^SIN  [AL]  ]  «-CLALT«ST«C9S  [At,]  /AREA 
CLTB  =  SIN  CaL] *XKWbT*CLALT«ST^C0S [ALI/AREA 
clbt=cut-cltb 

CLTD9  =  XkTB*CLALT*S IN  [DELTA]  ^ST^C0S  [AU'i-DELTA]  /AREA 

cltb*  [xktb+xkbt]  *clal.t^sin  [Delta]  ^^st«c0S[al+deeta]/area 

CL3DT=CLTd-CLTC3 

CLVISt*  C  [sin  Cal  +  PELTA] ^SINTAL  +  DeLTA] I »ST*CeS [AU  +  DeLTA] /area] «PDCT 
CLT'ClT+ClVIST 
6095  CC9T  =  c'^P*[ST  +  ExS]/aREA 
STTST=ST+EXS 
XL AYT4=XLaM14 

IF  CIz7,Y"4]  1610/  1610/6098 

6  098  XCPTB  =  XtAiL  +  C  CxKwBt*SIN [AL] ^XBCrWB  +  xKtB#S I N  CDELT A] »XBCRWB] / 

1  rxK^NBT-^SlN  CAL]  4X)<tb«SIn  CDFLTA]  ]  ]  »CR00TT 
8C  IF-[]2ZY  ^  4]  l6lO/  1610/  1710 
1C  XL09  =  X-l/D 

2xk*Vx^*ABS[SInCaL]] 

IF  [Zx""-'«8i  1310/  1350/  1350 
1-lC  CCC®2.4»SQRTCl,5i2g»1.5l29«ZXN®ZXM] 

98  T9  1391 
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135C  IF  tZx^-1.15]  1380j 1370> 1370 
138C  CCC»1 .64.SGRT  C.c34i^- C2X^'-.975]  ••2] 

33  T0  1391 

137C  IF  !:Zx‘-'-3.]  1360,  1381/  1381 
136C  CDC«1.9-Sg9TC.36i-.09«CZXP-3.]**23 
30  TP  1391 
1381  CCC«1.3 

1391  ittAs  CC  •  0000^35  *  CxL  OB  *«31 3  -  tO.  C0i73»  cXLGb#*?]  ]  +  CO  •  0298«XL0B3  +0.514  6 

TF  1395,  1395,139? 

1392  IF  CVX^^-I.A]  1399,  1394,1394 

139  3  ET^'*FTA+  Ci."ETA]*CVXI^-.53#1  .lUllll  mill 
30  TP  13=5 

.  f  r-  ^  A  ^  i 

^  S'  i  -  A  ♦ 

1395  IF  CXl"B-1C.]  i3?o, 1330, 1340 
13pc  xt<2<l  ®-0. 0054  »  cXLPi3-»*2]  +0.1  04 ^XlRB  +  C * ^37 
GO  TP  16CC 
133C  yK2Kl=0.939 
■30  T P  1  600 

l34C  =0. 93-9+ CO.  0C1525*  [XEen- 10*  01  ] 

'160C  alP'AL 

IFtAt]  I6c2, 16cl, 1601 

160?  CrX'="C2C 

1601  [xK2Ki*SIN  c2.*ALP3  *C9S  [  ALF/2  •  1 1  »3 , 14  159*0^0/ [4  , +ARFA] 

l  +  ETA  *crC+[  CXL*r;]  /  AREA]  +CCSIMCALP1]  «423 
VC  =  ><OG/0 

C!^31«  [Xk2k1«XQ»Sin  [2.*ALP]  +CeS  CaLP/2 •  1 1  +3 •  1 4 l59»p«D»D/ 

1  14  «  +  ^ARPA»xPETl  , 

[PTA  +  CCC*  C  I  XL +  03  /  ARE  AT  +  t  [XcG-  txU/2.]  3  /D3  ■«■  C  CSI\  CALP3  **23  3  3 

l.LVXRE'" 

C^B®C05l+C^23 
l6lC  R[:=HrFT*Vx^^*XL 

IF  IPE'1«E063  l6ll, 1612, 161 2 

1611  4A®*C835 
v\3. «  «  .  2  1  1 
33  TP  1617 

1612  ?F  C'^E-1'E073  i6  1 3, 1  6 1 4, 1 6 1 4 

1613  aa=.C52 

V  */■  =  ".  177 
GO  T«  1617 

1614  if  [PE'I’EOB]  1615,1616,1616 

1615  AAS.C33 

y  \  \  .  1  4  8  g 

G0  TO  1617 

1616  AA=*C221 
yNN»-.127 

1617  r  Fn0r  =  A  A*RF+ »X'.N 

CD9B  =1  .G2«CFB'-X*  Cl  .  +  1  ,5/ [XLOB**!  »5]  +7./  tXL8B**3»3  3  «SSU3S/AREA 
IF  CIZZY-43  161S, 1618, 1619 

1618  1  .  1999999 

IF  C3w3  1620^1620,1621 
16?C  '-r;P'A  =  o* 

1621  IF  ISwP]  1622, i6p2, 1623 

16  2  2  ‘f  .’'P'-f'  2  s  C « 

1623  IF  CSt3  1624,1624,1625 

1624  rCOT^O* 

1625  i-eNST=CD8'R  +  C05A2  +  CDeT  +  CD03 
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ne  TP  1639 

u-is  ir  CVX"-1.2I  1631»16.31#163;> 

1631  IF  CVx^^-tg]  1638/  1638/1639 
1638  CDALZsCCea 

CC3''*''2t  =  CDgB+CDpW+cC6T+CO0W? 

30  TP  1663 

163S  IF  rs''0  1700/  1700,  1701 
17CC  CCDOS'^''  =  0.0 
S.'JTPtsC.C 
CCOW'  =  .0,0 
oa  TO  1702 

17C1  XX^'"  .  VXm*S3Rt  [cos  [XLAHW41  ] 

ccf''lTc*SGFT  CAas  [  CXXH*XXr"l  -1  .03  3  /  CT9vCW**0»  333333 

ATC'AR'i'i#  C  TO  VC  W«*,  3333333 

TZT'l 

164C  ip  [AtC-1.53  16^1/1642,1643  ■ 

1642  IF  CVx^'-l.]  1643,  1643/  1644 

164  3  FlXCT  s3.  3o81.-1.8S779»SO.^.lTr+ll*0916*SOMl7G«SQMTTC-  l8.60S7* 
lSG'^nc»«'3H.7,4633*SQMlTC**4  '  '  . 

CO  TO  1650 
1644flXCT=3.4 
GO  TO  165c 

1641  IF  CATC“.53  1645/1645/1646 
1646  IF  [VX‘'''-1  .3  1647/  1647/  1648 

164  7  FL,XCT  =  ?.4  79l7-1.42798^tSa.i^ITC.-.3244G5XSQMlTC*SQ'‘'ITC 

ne  T-O  165c  • 

1648  FC^CT B 2 . 5«.ATC 
G8  TO  165c 

1645  IF  CVx.”''”!.]  1649/  1649/  1651 

1645  FUXCT.  .5325»  .47202»SGMITC-c,0863i*S.0NlTC*SQ^‘ITC 
GO  TP  165c 

1651  FUNCTs. 55917+. 33'?33*SQ^'ITC-.671429*SCMITC«SQMITC 

1650  IF  CIZT-2]  1652/1653/1654 

165  2  nCDOSw  =  F  UMCT«  CT0VC'/-**1 .666673  *  [[GBS  [XL  AMW43  3  **2. 53 

rcc‘3w  =  pcre5w 

17C2  IF  Cs'.'i23  1703,1703/  1704 
17C5  DCCOS?  ®  0.0 
S’vHTqt  =  o.C 

CDeW2s  0,0 

GO  Ta  I7c5 

1704  IZT  =  ? 

yyy  sVx'^^SqRT  [CaS  CXUAa'243  3 

SCXITcsSCrT  cabs  [  [XX6'*XXM3  -1  O  3  /  fTeveW2*«0  •  333333 
atC'ARWj* crevcua**. 333333  ''  ■  '  ■ 

G9  TO  ,164c 

165  3  FC00S2*FU.nCT*  C t9VCx2* »  1 . 6666734  ECCas[XLAX243  3  *  +  2,53 
FCO0X'2  =  CCoaS2  ■  - 

1705  IF  CST  ]  1706,1706/17^7 
17C6  CCDOsT  =  0.0 

sttbt-c.o 

CDBT  »  C.O 
G9  T0  i7c8 
17C7  JZT  *  3 

Xx5«Vx"'*SqRT  CCaS  [XLAf^T43  3 

SC4'‘ITcsSGrT[A0sC[XX4'*XXK3-1  « 3  3  V  CTevcT<»»Oi>  333333 

ATC'Art* Ct6VCT»+, 333333  . 

G9  to  164o 


1 65^ 
I7c£ 

1711 

17CS 


1661 

166c 

166F; 

165E 

1666 

166? 

1667 

1664 

165? 


1655 


1632 


166  3 


DCD0St»FUnCT«  Ct9vCT*«1.66667J#  [tees  tXl.AMT43  3«*2.5] 
ncDST-DCD0ST  ^ 

CeVC  »  1 ,31  a  1 3-0  •36633*VX^^+•06638•vXM»*2*•  00601  *VXM**3+»  000275^ 

CDFPTrxi #02«CFE9D*C8VC»SSUSS/AReA 
CDPPTR»CDeB«l .O3*CFE0D»SSuBS/AREA 

IFCVXM-1*C]  17C9> 1709^ 1711 
CCP2TR* cCDPPTR/0.21 » Cl .2-VX-3 

rR=XLNQS2/D 

CDWNlaO.O0O4O7>CFR»«8]  -O*  01  C2»  frR**?^  •i-0.l08*CFR*-»61-0.616*[FR«*bJ 
1  +  2,074»  CFr»*4]  «4. 183»  C^R*^^33+4.891^^  [FR*«23  *3 1 017  *FR  +  0 . 7795 

CC'*''N2aO.0c0l72*CFR«*S3-0»O0453*CFR«-»7]+0f050^*CFR**6]'-0«3C4#CFR-»«53 

l+liC96*CFR**4]»2t4C6*CFR*#33+3»l60»cFR*«?3"2*391»CFRl+1.0CO 

CC‘^^'3*.000l25»  [FR  **  83  «.  00370#  CFr*#73  +  ♦  0447*  CFR**63  -  .288*  CFR**53 
1  +  1.076*  CFR**43  .2.385*  CPR**33  +3* l41«  CFR**23 - 2 . 529*FR  + 1 . 30C 
IFCVXM“C*81  1664.1661^1662 
CCP7R«CDWsj1 
G0  T0  1658 

IF[VXm-1.C3  1655.  1665/  1666 

CCPTRsCDwnP 

G0  70  1658 

CcPTRs  C  [CdWN2-cDWN13  /0*23  *  rVXM'’0*83  +CDWN1 
GO  TO  1658 

IFCVXM°1.23  1667/1668/1664 

CDPTRxCDWn3 

GO  TO  1658 

CoPTRx  C  CCdWn3 -cDWN23 /O .23  «  rVXM-i .03  +CDWN2 

GO  TO  1658 

CCPTReO.C 

CCeBBC0FPTR4CDPPTR+CDPTR*VVH«3»l4159*D«D/C4»*AREA3 


CCAUZxCD9B 

CD9Ts  1 ,  1*  [DCOOST*  CSTTOT/ARfA)  -^CDOT] 

CDOwal , 1* [DCDOcW*  CSWT0T/ARrA3  +CO0W3 
CD9^2sl ,  1*  CCD0w2  +  BcCes2*  CS'-'TOT/aREA]  3 

CCO'-^HtM  .1*  CCDPW  +  DCCPsi-.*  rsvT0T/AREA3  +CD0W2+DCDOS2*  CSW2 TOT/ AREA 3 
l  +  CD0T  +  PcDeST*  CSTTQT/ARFA3  +rD0B3 
IF  IIZZY-A]  1659/1659/1663 
TONSTxCDOwBT 

t-0NSTsHsNsT  +  ZF*  CT0NST.HeNsT3 
yyMsVX^Rl 
IZZV  X  I^ZV  +  1 
GO  TO  1 

CD«^BT=  C  CH9MST.T0NST3 / CSQRT  C3.3 -SQRt  Cl »23  3  3  »SQRT  CVXM3  +TONSt 
1+  CM0KST-T0NST3 /C1..SGRT[3,3/SQRtC1.23  3 
COOT*! ,1* [DCDOt  *  [STT0T/ARFA3  +CdOT3 
CDOW*!  ,  1*  [DCDOiv  *  CSWTOT/ARfA]  +Cd0W3 
CCO'^^P*!  .  1*  CCOeu2+DCDOW2*r5'‘'T0T/AREA3  3 
CL'HBaCOOWBT-CDoT-CCeW-CDOW?'^ 

CDALZsCDOE 

CCb“XK2Kl*SIN  [2»*AL3  *S I N  t AL/2 . 3  #3. 1 4l 59*D*D/ [4.* ARE A 3  +ETA»cCC*XL 
1  .D*  C  rslN  m3  3  **3l  /  area 

CLB®Xk2k1*  CSIN  c2,*aL3  3  *C0srAL/2,3  *3, 14159<^D*D/  C4.*ARrA3  +ET4*CDC»XL 
l*C*t'^SlNCAL33**23  *ceSCAL3  /  ARE  A.iC0ALZ*E8S  CAL3  *COS  CAl3  *S  I  CAl] 
CA5*CD9wBt 

JFCAL3  1603/1604/1603 
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160'!* 

1603 

16C5 

1606 

1607 

160S 

1S71 


1S7C 


1S6S 


181C 

1820 

183C 


18^0 

185C 

18CC 


237C 


VCPBaO'C 
G9  T9  1605 

XCP3a  ctxCG/XPEr3- CC^'B/CNB]1#XREF 

IF  CSwl  1607/1607/1606 

TF  CSj]  2973/2973/1608 

IF  CSy]  1900^1990/2974 

F=D/2. 

IF  CICSC  -  n  l97C/l97Cil971 
xqi-p.w/b* 

XB2*BT/2* 

TJaD/BT 

MW1'-0‘5«cR00Tt*ABS  CSIK  CDfCtaI  3  +  tXCpWB»XTAIl.*CReeTT3«ABS  tSiN  CAt3  3 
XLAXlaXUAMW 

GO  TS  1969 
XEl “9T/2  • 

yB2*Pw/,2* 

TT=D/B'''' 

MWl*  CxCPTB-fXWlKG^CReOT'O  »  AE?S  tS  I  N  C  AL3  3 
X'L  Af^l  eXOAmT 

FTP7a[[xB2-R]/c2,#!;i.-TT3]1»r  13#  H 159/4.  ]  -  [  13*  14159#TT*'«23  /4*  1  •■TT  + 
1  [  [  Cl  .  +  TT**2]  /  12  ,«  Cl  ,-Tf'**23  ].]  *ArS  IN  1 C 1 . /  1 1 . +  TT»  »23  3  3 

FwsFTRT+R 

F  II  '  1F'>''«R**23  /  c  1fW#-»23  +  IMWI  *  «-23  3 
FI  1=  ChN1«R#*23  /  1  CPx«*23  +  CHu1*«-23  3 


20=hwi 

zlt=o.c  .  ■ 

00  1800  1*1/4 

ZLl®  1  IXBl-  iXLArUR]]  CZC^»  11  «-XLAf^l3  3  3  /C2,»  tXBl-.R3  3 
71,2«Al.9G  1 C  1ZD»»23  +  c  CZC-XBll  *«H3  3  /  C[zD*f»23  +  C  C2C-R3  *-*23  3  3 
ZL3®  C  Cl • "XLANii /  iXBl"R3  3  «  C  fXBl-Rl  +  CZO«  CAJAN  C  tZC*XBl3 /ZD3 -ATAN  C [ZC- 
1P3/ZR33]] 

ZL=  1^7L1*Zl23 -ZL3 
IFCI'21  181C/1E20/1820 
ZC“-ZC 
G0  TO  18.50 

IFCI«3]  lg30/ 1840/ 1850 
ZL'^'ZL 

ZC  =  2I1 
ZD^HI  1 
G6  TO  185c 


ZL  =  *-ZL 

zc--^ii 

.7LT  =  ZlT+Zu  -  . 

CONTI NUE 

IFClCSC-1]  297c>?97C/2971 
apt = CBT-DI *»2/GT 
XB  T  *BT/2 • 

r.Ll*  CCLALW«CUALT*XKWBl^SlNr  AL3  *SW*2.  *ZtT*  [XBT-R3  3/12 
irTRT«AREA»  C1»+xI-AMT3] 

CLl^CLUCeSlALl 

clt=clt+cli 

clit=cui 

rLlW=o. 

XCPTVeXCPTB 
GO  TO  2972 


.*3.14159*»aRT» 
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2971  ARW* CbW-DJ »«2/SW 
XBW’Bw/S* 

CLI*  faALw*CUALT*CxKWB!»siMCAU3+XKTBl«SlNJCDEUTA]]»ST«2.#ZLT«  CXBW- 
1R3  3  /  C2*»3, 1^155* ARW#FTRT*APEA*  i:i»+XLAMW3  ] 
rLl»CLl*C8SCAL3 

clw»clw+cli 

CLl^^CLI 

CLn»0. 

XCPTV»XCPWB 
G8  T3  2972 

2973  CUT'O, 

CLAUToO. 

CLTDbq. 

CDTD»0f 

CLBT=0» 

clbdt.o. 

ClTBbq. 

CLTDB-0. 

CL'^IST*C» 

CLI'T'O' 

CLlWsQ*  ■  . 

CLl“0, 

Ge  Te  2972 

2974  CL^'C. 

CLALWsO* 

CLB^'O* 

CL^B^O* 

CLVlSWsO* 

CL  1^*0* 

CLIT»0» 

CLi'C, 

2972  IF  CSw2]  700i7cC/7cl  .. 

70C  CLW2=0* 

CLALW2»C» 

CLB'^^EbC. 

CL^B^bO. 

CLVI W2»0» 

7C1  ALPPAbAL 

IF  CVXM«,0»53  1975,  1976/  1976 

1975  XK®1*17 

Ge  T8  1973 

1976  IFtVXM-l»3  1977/1978/1978 

1977  XK=2*0-SaRTC0.764-[VXH-0*1263 **23 
G0  T0  1973 

1978  IFCVXm-2*03  1979/1979,1980 

1979  -irK-^e^'SnRT  Cl  .293*  CVXH"2»13  3**23 
GO  TB  1973 

198C  yK«C»87 

1973  CDT®XK*ABsCSINCAL33*ST/AREA 
CCW'=XK*ABsCSINCAL3]  *SW/AREA 
CDW2  =  X'<*AbSCSInCaL3  3  *SW2/APEA 
CDD®Xk*a3sCSIN  CAL  +  DEUTA3  3*ST/AREA 
CCTObcDD*-CDT 

CATD“CDTD-*CeS  CaL]  ^CLTD*S  I  N  r  AL3 
CKTD  =  CLTD«C8S  CaL3  +CDTD*£ I N C AL3 
CKT»CLT*CeS  CALpHAl  ■*CDT*SINCALPHa3 
CAT^CoT^CeSCALPHAl «CLT*S I N C ALPHA3 
CAW2«CDW2«Ces  caU3  ®CLW2*S I N cAL3 
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VO 


CA'w"CDW*CeS  CalPwa]  •CLW#sisj  CaLPMAJ 
CNW»CLW#CeS  CAL^WA] +CDW#SrM  fALPHAl 

CN'a'2«cLw2«C8S  CaL] +CDW2»SlN  rAU3 
CAO»CdB#C0SCAL]-CLB*SI\CAL] 
rA5«CAT  +  CAW  +  CAl\2  +  CATD  +  CAB 
^0  T0  1991 
199C  CN^CN'B 

ycP2»xCPB 

CA=CA5 

CL’^Qf'CLB 

CDTST*CCB  +  CA»CPSCAl]  »C0SCAr3, 

CLALTsO. 

CLALWaCf 

roT'C. 

CCD'C.O 

CCW»G, 

CAT0=0' 

CMD  =  o. 

CKT»C, 

CAT»C  , 

CAW«C, 

CNW=C, 

CL-'J  =  <^. 

CLT'C, 

CLI*0» 

riTT^o* 

rLwTaQ* 

CLALW2»o* 

CPwBso* 
rAV.2  =  o* 

CNTT=o* 

CK'-'Tsc* 

:\KT2=C* 

CK'A'2aC* 

CL-v2  =  C. 

CL-‘-2T  =  o. 

XC’^T=.c* 

G8  T0  1992 

1991  CK=CNB-f-CK^(+CNT+C,VT04.CN]W2 

CLTT®cLET  +  CLT6  +  CLB0T  +  cLTDB*a  iT  +  CLVjST 

B8SCAL3-CLT0T*SlNrAL3 

IF  CST]  BBOl/ogoi^aqpP 

lVCl^T'''  +  rinI-,-^,y!:^p^^'*^I^^*Sl'^'CAL3]  »XCpTB+CL8T»C9s  CAL]  *XCPBT  +  CLTDB 

S21  IP  ^^LlT»C9srAL]*XCPTB-CLVlST*CasCAL].yCP?^^ 

i.b’'-cJ  9903,  9903/  1993 

7h 


S9C 


1993  ycPW2.  CtCLWB2*c8SCAL3+CDW2«SINtAL3  3  *><CPWBe^CLBW2*CeStALJ  *XCPBW2 
1  +  CLVI^<2*C0SCAL3«XCPWB23/CNWT2 
9SC3  IF  CSW2:  1994/ i9«54^  1995 
199^  VCP^2mO« 

CNWT2«0. 

TF  CSw3  1997/  1997/  199-5 

1995  ycP^=CCXCpW  ]  *CNWT+ CXCPWS  3  *CNK'T23  /  CCNWT*CMWT?3 

CN'-'JT  »CNWT  +  CNWTc 

IF  CSw3  1997/ 1997/ 1998 

1997  ycPWaQ* 

1998  IF  CST]  1999/1999/1996 

1999  ycP'*’’©* 

1996  ycP2»  CC^B*yCPB+XCPT*CNTT■t•XCPW*C^WT]/  tCNB  +  CNTT  +  CMWT] 

1992  CxCg2''XCP2] /yREF 

7CC1  AL1''*L*57, 29578 

87CC  print  5G02/  AL 1/ CLTei/ CDTet/ CL''*<T/ CLTT/ CL0/ CL  I  /  CNWT  /  CNT  T/ CNTC/ G\B 
lCN/CA^yCPW/XCPT/yCPB/XCP2/rP 
5CC2  FBRP'At  ClX/F4.c^  l7Ciy/F6.2l3 
ALPHAsyAL  [J+13 
6C0C  reNTiNUE 

print  94Cc^  CLaLT/CUALW/ClaLW2 

94cC  F0Rf^ATC///2Xj  6hClALT« /  r6  *  3  /  5X/ 6hCLAU^» ^  P6  •  3/ BX /  7WCL ALW2  =  » *  33 
''ELTA1  =  xDt  Cl  I +  l3 
6C01  cbntimje 

VXN'XvXP C lJ+13 
6CC2  CBNTInjUE 

IF[NPeDY-lL3  3334/3334/3333 
333a  CSNTIN'UE 
ST9P 
end 
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SLSR9UTINE  GE9SUB 

Clf'ENSieN  XVXMCl63/XDTtl6]  |XALt483 
XVx^l/XDT/XAL  ' 

re^MaN  CN,CAiCKB,C;'JW#C''iTiCvW2/CL.W/CLT#CLW2/CLl/CLW8iCLVISWi.CL.lT^ 

.  1  CLlW,CNTD,CAB,XCPaiyCG2>xCPB,XCPTiXCPW/XCPW2/XCG/XCPTV> 

2  XLa'^TA 

Cn^^^^9N  LLKK,  ISl>.PW/  ISWPT/  ISV'PW2i  f  AFBW/  I  ATgli  I  APBW2^  IL/LLLL>  IJ/ J# 

1  I  I,N89DY^  IZZVi  ICSC^  T^J8SE,NhiNHU</ IDTy  IM/ lAL^  ISWPI/ I  AFB 

XUA’^WiXLAyT^XLA‘^W2^yMACW^XMACT/XMACW2^CUAMW,CLAMT/CUA^W2/ 

1  BW^BTj  Bu2/CRe9T''^'/CR99TT#CR99w2iSW,ST/SW2,XWlN5iXTAlLiXWlNG2 

2  ^XLiDiDi/XLN/AREAiXREF/SsUBS^XLSB^ZF/ARTiARWMRWB/XUNeSE 
CSM^'SN  CeLAM^BC9LAM^CR9eT^nl,CLAUliXl.AHl|BAR/RATI0,X'<TB/XK3Ti 

1  XKwB/XKpW,  XDCRBW/  XK'/BWjiXkWBT^  XKWBW2  /  XKBWW^  X<BWTj  XXBWW2# 

2  XCPBWiXcPBT,XCPBW2^vCPWB,XCPTBiXCPWB2/XKWRI ,8DC/9DCW/9DCT , 

3  eOCWEi  ClALW,  CLALT  /  CI.At,W2,RE'FT^  BETA  /  At  /  T9NST/H9NST/HT  /  XKTBI/ 

4  VXM/VXMr1,DELTA/XKCPWB/XbCRWB#XMAC^ XLAMW4/XLAM24/XLAM2/ 

5  XUaM^4  TeVC/TeVCT/T9VCW2/TBVCWiEXS/ STT9T/ SWT9T/ SW2TST/ RE/ 

6  CD0^CD8W^CD9T/Cn0B/CD9W2/CLBW/CLWB2/CLBW2/CLTQ/CLBT/CLTD* 

7  CUTDB/CL3DT/CLVIST/CLVIW2/CMB/CDewBT 
C8^'M9^;  DCD9SW/nCn0ST/DrO9S?/DCD0W/DcB0T/DCO8W2 
Ce^'^"9M  CDaUZ 

CRAT  e  XtNeSE/O 
ni”9 

Xt0B  *  Xt/D 
IF  rCLA’^'A"50  IC^  10/20 
1C  CLAf^W  =  5. 

2C  IF CCLA"w2-5.3  30/30/40 
30  CLAf">VJ2®5« 

40  IFCCLaXT”50  50^50/60 
50  CLAf'TsB, 

60  CQf^'TlKUE 

CLAf^WsCLAMW/B?,  29578 
CLAHT=CLAmT/57. 29578 
QI_A^<W2sCLA^"W2/57.^9578 
IF  CINSSE^U  1571/  1572/  1571 

1572  SS1  =  8.*3.»  [DRAt-2/51 
nn  TP  1579 

1571  IF  ClNi9SE-21  1573/  1574/  1573 
1574  SS1  =  4.4.f,*  CDRAT'’1«5]/3* 

^0  TP  1 b79 

1573  IFClNeSE-33  1579/1576/1579 
1576  SS1®4,+2«* CDRAt-I/S?] 

1579  sSUes»SSl*  C3.  l4l59*0»D74,]  +3. 14i59*D*  [XL-.XLN9SE] 

IF  CDRAT-.6]  1580/1580/1581 

1580  ZF=1* 

00  T9  1594 

1581  zr=l*522®EXPC-.7»DRAT] 

1594  z=HT*0*OCi 

Ijr  [HT-35332.0]  1595/  1595/  1596 

1595  Te5l9,.>HT/280. 

pS» Cl*51-0*01315*Z] **5»256 
00  T9  1597 

1596  T»393, 

B® 1 ‘SS-O* C478«Z 
PS=6»49»ExP  C93 

T6 


1597  C*'*9»1»SQRTCT3 
PS®PS#70»  9 
RH8"Ps/ CI7I5  9«t3 
XNU^2.270*ET»*1s5]/C  CT<'198»63  *fiO.«*83  3 

rept» [c*RHe3/XMU 

return 

f\d 
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SUBR'^UTIN£  CLASU9 

PI^'ENSleK  >;VXMtU]iXDTi:i63,)(ALCA83 
XVxMi  XDT/ XAU 

CQMMeN  c'^iCAi  CK^BiCNW/C^'T,c^W2,CLW/CLT^  CLwa#CUlj  CUWB/CLVISW|CIIT/ 

1  CLlW/CNTD,CAB.XCP2>yCG2iXCPB,XCPT^XCPW#XCPWa.XCG/XCPTV, 

2  XLaMTA  ■ 

LLKXi  ISwPW#  ISWPT*  IS^‘.PW2>  IAFBWi  lAFBT#  IAFBW2,  IL/LLLL/  IJfJi 
1  ■  I  I^NBBDy#  IZZY/  ICSC/TN8SE,NM,Nf^lX#lDT/  IALi  ISWPl  W  AFB 

XlA'^lW  iXL  AMT/  XLA-'*W2^  y^^ACW>XMACT/XMACW2/CUAMW|CLAHT/CLAMW2, 

1  BW,BT/Bv-2»CW8eTW#CRPeTTiCP80w2#SWiST>SW2/XWlNG#XTAlL/XWING2 

2  ^XL/D4DliXl.^J>  ARFA^XREF/SsUBS/XUSBiZF/ ART/ ARW/ ARW2/XLN9SF. 
rgMYoK.  cPlAY/ BCBL AMj  CRP0T^  n 1/ CLALI/XEAPI/ BAR/ RAT  18/ XKTB/ XKBT t 

1  "  XKw3/XKBW/XBCPB'^./XkvBW/XKWBT/XKWBW2/XKBWW/XKBWT/XKBWW2/ 

2  '  XCpEW/XcPBT/XCPPW2yXCPWB,XCPTB>XCPWE2/XKWBI/8DC/SDCW/0DCT/ 

3  PDCW2/Cl.AlW/CLAI.  T/Cl.ALW2/REFT^BFTA/AL/TeNiST/M9NST/HT/XKTBI/ 

4  VXR/ VXMRl,  delta /XKcrWB/XbCRWB>  XM  AC#  XLAMW4,  XLAI^24/ XUAM2/  - 

5  XLaM4/ TpVC/TBVCT/TevCWE*  JBVCW# EXS/ STT0T/ SWT9T/ SW2T8T/ RE / 

6  CD0/CDB‘A/CD0T/CD8B/CDBW2/CUBW#CLWP2/CLBW2/CLTB/CLBT/CLTD/ 

7  CLTDB/CLBDT/CLV!ST^rLVlW2/CXB#CDewBT 
reYPBN  CCDSSW/CCDesT/DCDPSP/DCOeW/DcCeT/DCDDWB 

cdalz 

IFt:VX,M-l.]  2/2/3 
c  EETA=SSPT  C1.-vXY*«-23 

G9  TB  4 

3  peTA®S2vRT  CVXm**2-1,3 

4  lFnZZV-4]9C2/90?/llll 
nil  IF[VXM-1. 3  41/42/41 

42  PETA=o«CCo0001 
41  XFIN'C 

KFIN»kFIN+1 
IF  CS'A3  5C4/5CA/411 
411  ARw  »  [yw  -  0]*-*>2/SW  .. 

«AR»BeTA^*aPW  '  (  n  n  /■ 

iswPnsuPu  ,  i. 

XLAY=XLA^W 

ARsARw 

5C5  IF  ClSWP-n  5/F/20C 
t  IF  CXLax-,25]  AO/lO/10 
1C  IFCVXM"1*03  20/20/30 
20  CLAP  =  - • 1 S33*BAF+1  *6 
OB  TO  370 

30  IFCBAR-^1.]  50/40/40 
40  CLAR=1 .50s* Cl  .26] ** [2.-BAR3 
G6  TB  370 

BO  c L aR ~  •  3 ■» B aB"*"  1  •  6 
08  TO  370 

6C  TFCVXY”!.]  70/70/80 
70  CLAR  =  ”*l667*3AR-f'l.575 
Gfi  TO  370 

8C  IFCBar-1.3  90/90/100 
90  CLAP*. 1667*BAR+1 .575 
OB  TO  370 

ICC  IpC3AR-2*3  llO/llO/lPO 
lie  CLAR^I  •74i7 
OB  TB  370 
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12C  *^28*  tl  .223  t2.-BAR3 

G0  TS  370 

20C  IF  CXlAM^.I]  2l0/260#2fc0 
2K  IFCVXm-1.3  220*220/230 
22C  CLAR'!-»2077»BAR+1.575 
Ce  T9  370 

230  IFCBAr-,253  24C<240/25n 
240  CLAR»  .2o77*BAR  +  1 .575 
G8  Tfl  370 

250  IFCBAr-4,]  251,251/252 

251  CLAf^*-.*  1668*BAR+1 .667 
G0  T0  370 

252  CLAf^  =  1.587»C1.26]##C2,,BAR3 
G9  TP  370 

260  IF  CXLAM-,3]  270/320/320 
27C  IFCVXk-I']  280/280/290 
280  CLAR'-*2065ftBAR-*'l  .6 
G0  T0  370 

cSO  IFCBAR-C.753  3C0/300/310 
300  CUAR'.2o65«BAR+1,6 
G0  T6  370 

310  If:BAr-2.53  311/311/312 

311  CL^^*”*217«BAR+2.293 
G0  T8  370 

312  CLAB  =  1 .543* Cl .263 **  [2.-BAR3 
G9  T0  370 

320  IFCVXm-1*]  330/330/340 
330  CLA^»- •?25*BAR+1 .675 
00  TP  370 

340  IFCBAr-1.]  350/350/360 
350  CLA«».225*BAR+1.675 
GP  TP  370 

360  CLAR*l  .508*  Cl  .26]  *«  C2.«.BAR3 
37c  ifCAR-I.C]  800/800/810 
800  ARR-I.O 
GP  TP  820 

810  aRR=  1 ,/ C-Ar**  C  CAR..1 .3 /AR3  3 

S2C  CLALscLARttARR#/iR 

■IF  CkFIN’.2]  500/501/502 
5CC  CLALWbCLAl 
GO  TP  503 

501  CLALw2=ClAL 
GP  TP  503 

502  CLALT=CLAl 

5C3  IF  CKFlf."2]  504^507/900 
504  ^FlNa|<r  yN+1 

IF  CSw23  507/5c7/506 
506  aR= CSW2"D] **H/sW2 
3AR'“BeTA*aR 

ISWP=lSWRw2  5f'f 'if'  Ql  do  ^ 

ARW2saR 
XUA^^•‘XLAMW2 
G0  T0  5c5 
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5C7  KFI^J'-KFIN^.I 

IF  tST3  5o9^509^508 
5C8  AR« CBt^DI #*2/St 
BAR*2ETA*AR 

ART*AR 

iswp«iswpt 
VLAMaxLA^‘T 
G0  T0  505 


5C9 

IF  CSwl 

5100^  5100/ 900 

51CC 

IF  CSw23 

5110/5110/900 

511C 

IF  CSt3 

980/980/900 

98C 

LLLU«2 

return 

5CC 

llKK^Q 

5C2 

IZZY  a 

IZZY  +  1 

ir  tSw]  5l0^51o.»9Cl 
9C1  CSLA^'a  COSCCl,A^'W3/[SINCCL.AMW3] 

ARW*  tBl'i-D]  «*2/sW 
eC3l-AM»BETA*C9l,AM 
CReBT.CRPaTW 
B1=BW 

IAFBoIAFBW 

XRACaX'^ACW 

T3VC*t0VCw 

CLAUbCLAL'*-' 

yLAMl.XLAMW 

ISW^l=ISWpW 

par*beta<^aRw 

RATlP-CReeT/CBETA^OI 

LLl-L«0 

return 

51C  IZZY  a  IZZY  +  1 

IF  CSW23  9A2;94p,5H 
342  LLL1-  =  1 

return 

511  CeUA^'  =  C0S  tCUAMwB]  /  [S I N  CCi,  A»^W23  3 

nc9LAMaBETA*C0LAM 

CRQ9TsCRee’-^2 

PlsBVvE 

lAFBalAFBwB 

CLAUIsCLAlWS 

VLA^I aXEAMWB 

y V ACaX“ACw2 

T6VC*Tv‘'VCw2 

I  S‘'51 '  iSWpWB 

ARW2= [BW2»D] »*2/SW2 

par®bet  A^ARI'^B 

rati BaCReoT/ C8eTA«D3 

LLLlaC  1 

return  ‘  ' 

END 
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subroutine  catsub 

niNENSlON  XVXMCl63|XDTri63 iXALC^S] 


CfiNM.BN  XVxMiXDT^XAL  ,, 

CN,CAiCNB,CNW,CN'T#CvW2#CLW/CLTiCtW2iCLl  iCLWB/CLVlSW,CLIT^ 

CLlW/CNTD<CABiXrP2,yCG2/X,CPB/XCPT#XCPW/  XCPWSiXCG/XCPTV/ 

XLAMT4  ,  ,  ,  ,  , 

CBHNBk  LLkX/  IS'.',PW/  ISWPT/  ISI*'PU2/  IAPBW/  I  AFBT^  IAFBW2/  IL/LLLU^  lU^  J/ 

I  I ,NBBDYi IZZYi ICSCi iN6SE/NM,NMLK/ IDT/ IM> I AL^ iSWPi/ I AFB 

CeNNBxi  XUaNW/  XlANT/  XL  A''W2/yNACW,XHACT /  XNaCW?/  CUANW/  CLANJ/  CLAHW2/ 

EW^  BT/BK2,CR8eTN’/CRBeTT/CP0Bw2/SW^ST/SW2/  XW  I  NG/  XTA  I  L/ XW  I  NG2 
/XL/D/DI/XUN  /  ARPA/XPEF/SsUBS/XURB/  ZF/ ART/ ARa'/  ARW2/XLNRSE 
CQNNRNj  CRLAN/BreLAM^fRneT/Pl/CLAUl/XLANl/BAR/RATie/XKTB/XKBT/ 

XKw5/XKeA/XBCRBA/XKyaW/XK'A'BT/XKWBW2/XKBLW/XXBWT/XKB.-iW2/ 
XCpBWi  XcRBT/  XCP'^wa/ vCPWB/  XCPtB  /  XCPWBE  /  XKWBI  /  ODC#  9DCW/  8DCT/ 

eDC^'2/  Ci  ALW/  CLALT/  CLAI.W2/REPT/BETA/  ALi  TRMST  /HRNST/HT/XKTBI/ 

VXM/VXMr1/DFLTA/XKCRW3/XbCRWB/XNAC/XLAMW4/XLAH24/XLaM2/ 

XLaNA/  TeVC  /  T9VCT/  TevCW2/  TRVC'a/  EXSi  STTSTjSWTRT/SWETQT/RE/ 
CD0/CD9i.,/CD0T/CD0R/rD0W2/CLBw^  CLWB2  /  CUBW2/ CLT3/  CLBT  /  CLTD/ 
CLTDB/CL3^T/CL.V!ST/^LVIW2/C'^B^  CDBWBT 
CPNNa^j  CCDeSW/nCDesTiOCDGSP/DCn0W/DcD9T/DCDeW2 

CRNNRn  CDaUZ 

XKXB»C2./3.1A159]^^|;C1,+D«♦A/B1*^^'^3^^[»5^^ATAKC«5•■[B1/D"D/B13]+3.1A15 

19/A.]  -  Cc<^*?/B1«*23  C  CBl/D^P/ei3  +2**aTAN  CD/81]  3  3  /  Cl  .-D/Bl3  **2 

926  vKTBl*C[3.1A1593**2*D*»2-*CRl/D  +  l*3*-)t2]/[4»*3l**23  +  C3.1"’l5^*D*»2*CB 

11.^»2/B**2+1  .3  «*2/  CB  1^**2*  CGI /D-i  *3  **2^  3  *ArSIN  C  C31  *»2/D»  *?»  1 . 3  /  C51<^« 
2P/D**2-<-l'33 

XKie2s2.  «3. 1A159*d*  CBl/C  +  r*3  /  [Bl*  ^B1/D"1.33  »  C  CRl  »*2/D**2  +  l  *3  **2/  C  C 

lBl*»?/D^(♦2]*CBl/D-1.3**233»CARSINCCBl**2/D**2-1.3/CBl*«2/D^^«2+1.3J 


1 

2 

( 

1 

( 

1 

2 

( 

1 

2 

3 

A 

B 

6 

7 


23  *  t?  r  ■  o 

XKT33cCA.*D*  CO  l/D  +  1  . 3  /  CBl*  CBl/D-l  *  3  3  3  ♦ARslK  [CBl.i*2/D»*2-3‘3/LBl**2 
l/D**2+1.3]*CS./CBl/D-1.3*«?3«ALeGC  C3 1  ^*2/D^<  *2+ 1 , 3  /  [2,  *D  1/03  3 
vkTB=  Ci  C3.  lAiBg]  **2]  *  CXktBI -X(<TB2-XKTB33 
IF  ClcSC-13  970^970/971 
S7C  IF  CLLXK3  972/572/973 
972  XKXBI=C/8i+1 . 

(5^KTRi=xXt3 
V  TR  973 

9>l-'iP  CLLXK-ll  973/973/972 

97  3  nAREF  =  BAR*.Cl  /  +  XLAM1]  *  [  C 1 . /  CBET A ^CRL A N]  3+1.3 
IF  CPAPpP-A.l  r2A/928/°29 
925  xkBa»  [1 , +D/B13  ••*?»XK'aB 
GO  TR  ICOO 


9291/9292/9291 
'1.3  932/931/931 


929  IF  ClAFBl 

=  291  IF  CGCRLAV’-I.]  732/931/931 
931  vkB'XU  [ECpLAN/  [1  .ABCRLAyi  3  «  C  C  C  CfiCRL A^  +  1  *3  ^  C 1  ./RATI  R3  +BC0LAM3  /BCRLA 
U']  «■«?]  aRc8S  [  Cl .  +  T'l.  +t^CRLA'''3  «tl^/RATl  63  3  /  CBCRLAy+  CBCRL  A"' + 1 . 3  *  C 1 «  / 


CBCRLA 

2PATIR333  ,  .  , 

XKB'>'2=  CSQrT  CDCeLAP*  +  2-l .]  /CBCetAN+l  ,3  3  *  CSGRT  Cl  .+2. 

1 1  J  »  CSGRT  [BCaLAN-»«2-l  .3 /3C^LAP3 
22FTA*D3  3+S3RTC  Cl  .+CR0R'^/ 


ftr’Ti  t  j  j  I  uj.  »Tfi»*BETA  +  D/  CRR0T] 

♦  CBETA*D/cR08T3  **2*ALRG  C  Cl .+CRRBT/ C 

_•  »  m  m,.  f-  ,  t  KA  J  f  A  .  ^  Cl  I 


3  Cl . /PcRLAh] 

XK3W«  [?  .^^BCOLAY/  C3. 1  a  159* sort  CBcR1-A,m*«2»  1 . 3  *  C1.+XLAM13  *  CBETA*D/ 
1CK89T]  .(  [Bl/D-1  .3  *  [hETA*CLAL13  ]3  *  CXKB^  1-»’XXBW23 
G8  TP  loOC 
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932  f  CSC0LAM  +  Ct  ,*gc©LAM3  «BETA®D/CRt§'n  /BC&UAMJ  ©el  1 5+ C  tBC0L AM*n  # 

l  +  eC8LAf^3  eBE:TA*o/CRPST]  /BCQlAM]  **^5^2* 

XKBW2.  C  [Cit+BCOlAM]  «BETA«D/'CReS?3 /BC^IAMJ  eeSeO  #  5e  CAU9G  Cl  • +SQRT  CBCe 
.1LA^^/  tBC8LAM+  Cl  ,^BCaLAMl  eBETAep/CRSI^T? 3 3  eALPG  C 1 . -SQRT CBC01AM/  tBC®l,A 
aM+  Cl  •■«'BC9lAM3  «BETA«D/CRSST3  3  3  3 

XKBw*  C16**  C0C9tAH/n  »  +  BC8Ca^13  3  13 » 1 A lB9e  Cl  t*Xi.AM13  •  CBETAeD/ 

ICReai]  »C0i/D*-lJe  t3ETA«CLALl3  3  3  »  tX<BWleXKBW23 
G0  TP  1000 

529c  IF  CBCOIAMp.1.]  9295^929^^9294 
929't  IFCRAtIP-1  *3  9293i  93 1 1  >  93 1 1 
5293  RAT I S»  1  •  . 

D»CRPeT/9ETA 

3311  XKBl’’  [1  «-*-CSLAMi?BETA«RATie]  *e2#ARCaS  C  CbCBlAM'^’RATI  83  /  Cl . -J-CBlAMeBETA 

\K3H“BC6LAP'**2'»RATie«*2*ARC8S  Cl  ./BCeLAP*!  fBCGLAMeRAT  1 8 **2*SQR'C  CBC9L 
1AK**3-1 *3  *ARSI\ Cl ./RAtI83 

yK33»SGRT  CBC8LAM*e2'-1.3*ALPQCRATl8-i>sQRT  CRATI0»«2*1.3  3 
yi^3>i=  [  C8./  C3. 14159*SGRT  CBCPLAr^**2- 1  .  3  •3ET  A«CL  Al,  1  *  CXL  AMI  + 1 » 3  #  CBl/Dl 
1-1 .] 3] *  Cl ./RATI  9  3  3  *  CXKB1-XKB2-XKB33 
GP  T9  1000 

9255  iFCRAjie-it]  9996/9312*9319 
9296  RATIPsI . 

rsCRPOT/OETA 

9312  XKal=  Cl  .-t-CPLAM^BETAeRATlR]  eSQRT  [  CRATIQ"!  *3  e  CC9L  Am*BETA*R  AT  I  9  + 1 . 3  3 

XKBc  =  RATI9«*2*C8C9LAM]  **  1 , R-BC^LAMeRAT I SeeBe  CBC0UAM4-1,1*CATANCSQRT  • 
l[l./Pc9LAM33'»ATANC3QRTCCRArI9"l,3/CC®l-AM*BETA«9ATI6  +  li3333 
XKB3aCCBC0LAM+i.l/SQRTCBCerAH]3eO»5eCAU9GCl»'i'S'3RTCBC9LAM*CRATI6*l» 

n  /CCeLAM^BETAepATlO+l.]]  3  -  AteGCl.-SQRT  CBCBl.AM«  CR  AT  I  8»  1 , 3  /  CCSLAM* 
2BETA*rAT  Ie  +  1  *3  ]3  ]  _  /r,. 

xkB^" C  C16.» sort  CBC9L ami  o  Cl ,/RATi93 / [BETAeCL AL 1 « CxL AM  1 +  1  *  3  *  CB 1 /D  1 
l-l.]*3.l4i59*CeCeLAM  +  i,]]33e  CxKb  1 ''XkBE«XKB33 
ICCC  yKBT=X<W3-XKTB 
OsDl 

IFCVXm-1*03  20CO/200C/2C1C 
2C0C  IFClSwPl-1  3  2020/2020*2030 
2C2C  IFCeAR-2.G3  2040/2050/2050 

20 4  C  XECRWe«=0*35-XLaMi«  [0 . 35+ 1  7 » 885"SQRT  t328  *  8782«>  CBAR-3  *  ]  «e23  3 
00  T9  ?C80 

2C5C  XBCRWB3C»35"0. 1*XLAM1 
GS  T0  2080 

2C3C  IFCBAr-2.0]  2060/2070/2070 
2C6  0  XCRl=SORT  C2029»5+  C3AF-3.3  e«23 

XBCRWb*  C  C-44.5  +  XCR1]  "XLAH'i*  4<t .  B  +  XC^l  17  *  885«SQRT  C328t8782-  CBAR 
1»3«]».»2]]] 

Ge  TP  2ilc 

2C7C  xaCRWB=C«55»0.3»XLAMl 
GP  T9^211o 

2C8C  IFCfcAR»4*C3  2090 » 2 1 00/ 2 1 OQ 

2 CSC  xpCRGw*0'25-XLA^l«  C32. 125-sQRT  Cl032. 02»  CbAR"4  •]  4*23  3  +  Cl  ."XL, AMU  * 
IalPG  Cl  .0^+0*  lep/BU  4  C-7.5+SQRT  C72»25«  C3AR"4.03  4  42U 
G9  TP  24CC 

21 CC  xbCRBw 4  0*25+  C  Cl .  -XL AMU  « ALPS C 1  •  O^-^O.  1«D/B13  3 
GP  TP  2400 
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2llC  IFCBAr-4#03  2l20j2130/?l3Q  ,  ,  _ 

212C  XBCRBw»0*25*  Cl . -2  .«Xl,AM13  «  (:32. 125-SqRT  tlC32  •  02-  [BAR -4 ,3  -**2]  3  + 
1  Cl  •''X|_Arl]  *AieG Cl  •  12  +  0.3*0/813*  C"7,5+SQRT  C72 .2B-  CBAR-4 .3  **23  3 
Gfl  Ta  2400 

213C  VBCRBw«0.25+Cl,-yLAM13*ALQnCl.l2*0.3*D/Bl3 


Ge  T0  2400 

2C1C  IFClSwPl-1  3  ?140.2140>2150 

214C  IFCBAR-3.o3  2160# 2l7Ci 217c  ^  , 

216  c  yBCRWe*!-  C9.235+25»«C1.-XLA‘'13  3+s0R3'  [  C9. 7 1 +25.  *  Cl . - 
1 CSAR-3.3 **23 
G9  TO  2200 

217C  X8CRWb»0.005*BAR*0.46 


ytAMl3  3  **2« 


09  TB  2200 

2150  IFCBAr-3.3  2180^2190/2190. 

218C  XBCRWb»0*675-XlAM1«  CO .  675+<» »  235-SQRT  C94.i-CBAR-3.3**23  3 
GS  T9  2300 

219C  XECR'^B  =  C.G05*BaR  +  0.46+0.2#  Cl  .-XLAM13 
GO  TO  2300 

220C  BARL,Am  =  BAR*  Cl .  -XL AMI  3  *  Cl  ♦  +  f  1  * /BcBL  Ari3  3 
IFCBArLAM-4.03  2210/2210/2220 
221C^rCBAR-2.c3  2230/2240/2240 
223cQ)x3RId1  *  AL9GC1.32  «  C.3e*XLAMi3 

2235  XB-lC2  =  4.  +  CXBMiDl«#23-CC0.5+0.5i39*  cD/B13  *C1»17  +  XLAM13*C1./C0.331 
1 ** CD/B13 3 3 3 **2] 

XB’-'IC3  =  2.*  C  CX3MID13  -  CO.  5+0.51 39*  CD/b13  *  Cl .  17  +  XUAM13  *  Cl.  /  CO,  331 
1 ** tD/Bl3 33 3 3 
XBMID»XBMiD2/XbMiD3 
yC^'I0»4,+  CXBMIc1-X8MID3**2 

XBCRPW=SGPT CXCMiD- CCBAR-2.3 **23 3+XBmID 

GO  TO  2400 

224c  XBCRBW"  0.5  +  0.256  95*  CD/B13  #cl*  / Co* 331**  CD/B  1333*8  AR*C1.17  +  XLA‘'' 13 
GO  TO  2400 

222c  IFCIAFB3  2250/2250/2260 
225C  IF [RATIO-I .03  2280/ 2280/ 2270 
228C  XBCRPW'0.67 
GO  TO  2400 

2270  VBCRBW  =  »2,32  +  SgRT  C8 . 94C 1 - C  C 1  * /R AT  1 9] - 1 . 3  **23 
GO  TO  24CC 

226C  XBCRBw* CCi429/RATl63 +0.5 
60  T0  2*00 

23CC  DARLAm=BAR* Cl .-XUAM13 * C1»+C1./8C0LAM^3 
IF CBArLAM,4.0]  2310/2310/2220 
231C  IFCBAr-2,c3  2320,2240/2240 
2320  XB'^101  =  AL6G  Cl  .65-0.65»XLAm{3 
GO  TO  2235 
240C  0*01 

arat*bar/bbta 
IFCARaT-3.0]  2C/20/21 

21  ARAT*3,o 

2C  hcC®2,*  C1,  +  XLAm13  »* 1 . 6*5X2 r- , 4* aR AT3 
IFCARaT«1,C]  22/23/23 
23  0DC®0DC-O.25* CaRAT-1 .03 

22  continue 

return 

FND 
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TABLE  2 


INPUT  NOMENCLATURE* 

AREA  arbitrary  reference  area, 

BT  total  tip-to-tip  tail  span  including  the  missile  body** 

BW  total  tip-to-tip  wing  span  including  the  missile  body 

,  total  tip-to-tip  wing  two  span  including  the  missile  body 

CLAIiT  leading  edge  sweep  angle-  of  the  tail*** 

CLAMW  leading  edge  sweep  angle  of  the  wing 

CLAI-rN2  leading  edge  sweep  angle  of  wing  two 

CROOTT  tail  root  chord 

CROOTW  wing  root  chord- 

CROO'R'G  wing  two  root  chord 

D  body  diameter 

HT  altitude  in  feet 

lAFBT  afterbody  constant  for  the  tail 

0  -  no  afterDooy  following  the  tail 
1  -  afterbody  following  the  tail 

lAFBW  afterbody  constant  for  the  wing 

0  -  no  afterbody  following  the  wing 
1  -  afterbody  following  the  wing 

IAFEW2  afterbody  constant  for  wing  two 

0  -  no  afterbody  following  wing  two 
1  -  afterbody  following  wing  two 

lAL  nvmiber  of  angles  of  attack 

ICSC  control  surface  constant 

1  -  tail  control 

2  -  wing  control 

3  “  canard  control 

IDT  nximber  of  control  surface  deflection  angles  (must  beat  least 

one;  if  there  is  no  control  surface,  IDT  =  1,  DELTA  =  OoO) 

IM  number  of  Mach  numbers 


INOSE 

ISWPT 

ISWPW 

ISWPW2 

I'JBODY 

ST 

SW 

SW2 

TOVCT 

TOVOT 

TOVCW2 

XAL 

XCG 

XDT 

XL 

XLMT 

XLAiMW 

XL.4J.IW2 

XLNOSE 


TABLE  2 
(continued) 

nose  constant 

1  -  blunted  ogive  or  cone 

2  -  pointed  ogive 

3  “  pointed  cone 

sweep  constant  of  tail 

1  -  unsvept  mid-chord 

2  -  unswept  trailing  edge 

swdep  constant  of  wing 

1  -  unswept  mid-chord 

2  -  unswept  trailing  edge 

sweep  constant  of  wing  two 

1  -  unswept  mid-chord 

2  -  un swept  trailing  edge 

number  of  configurations  being  run  (a  configuration  is  one 
complete  data  deck) 

exposed  planform  area  of  one  pair  of  tail  panels 
exposed  planfonti  area  of  one  pair  of  wing  panels 
exposed  planform  area  of  one  pair  of  wing  two  panels 
thickness-to-chord  ratio  of  the  tail 
thickness-to-chord  ratio  of  the  wing 
thickness-to-chord  ratio  of  wing  two 
missile  angle  of  attack  (degrees) 

missile  center  of  gravity  location  as  measured  from  the  nose 
control  surface  deflection  angle  (degrees) 
missile  length 

tip-to-root-chord  ratio  of  the  tail 
tip-to-root-chord  ratio  of  the  wing 
tip-to-root-chord  ratio  of  wing  two 
length  of  the  nose 
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TABLE  2 
(continued) 


XMACT 

XIvIACA^ 

XI4ACW2 

XREF 

XTAIL 

XVXM 

XWING 

XWI1'JG2 


mean  geometric  chord  of  the  tail 
mean  geometric  chord  of  the  wing 
mean  geometric  chord  of  wing  two 
arbitrary  reference  length 

distance  from  the  nose  to  the  leading  edge  of  the  tail  root 
chord 

missile  flight  Mach  nianber 

distance  from  the  nose  to  the  leading  edge  of  the  wing  root 
chord 

distance  from  the  nose  to  the  leading  edge  of  wing  two  root 
chord 


*  The  control  surface  is  defined  as  the  tail  regar^ess  of  the  mode  of 
control,  and  the  fixed  surface(s)  is  (are)  always  defined  as  the 
wing(s)=  See  Figure  lo 

All  linear  dimensions  are  in  feeto 


***A11  angular  dimensions  are  in  degrees 


TABLE  3 


PROGRAM  INPUT  FORMAT 


1 

Run  Identification 

Configuration  Title  and/or  Number  (1st  60  coliaims) 

2 

Control 

Constants 

ICSC  INOSE  IDT  IM  lAL  NBODY 

(15)  (15)  (15)  (15)  (15)  (15) 

3 

Wing 

Inputs 

ISWPW  lAFBW  XLAMW  CLAMW  BW  CROOTW  SW  Xl-IACW  XWING 

(15)  (15)  (F10.5)  (F10.5)  (F10.5)  (F10.5)  (F10.5)  (F10.5)  (F10.5) 

.. 

Wing  2 

Inputs 

ISWPW2  IAFBW2  XLAMW2  CLAMW2  BW2  CR00W2  SW2  XMACW2  XWING2 

(15)  (15)  (F10.5)  (F10.5)  (F10,5)  (F10.5)  (F10.5)  (F10.5)  (F10.5) 

5 

Tail 

Inputs 

ISWPWT  lAFBT  XLAMT  CLAMT  BT  '  CROOTT  ST  XMACT  XTAIL 

(15)  (15)  (F10.5)  (F10..5)  (F10.5)  (F10.5)  (F10.5)  (F10.5)  (F10.5) 

6 

Miscellaneous 

Data 

HT  D  XL  XLNOSE  XCG  AREA  XREF 

(F10.3)  (F10.3)  (F10.3)  (F10,3)  (F10.3)  (F10.3)  (F10.3) 

7 

Miscellaneous 

Data 

TOVCW  T0VCW2  TOVCT 

(F10.3)  (F10.3)  (F10.3) 

8 

Control  Surface 
Deflection  Angles 

XDT  Any  nmber  of  deflection  angles  up  to  l6  may  be  input. 

(F5.1) 

9 

Missile  Flight 

Mach  Numbers 

XVXM  Any  nianber  of  Mach  nmbers  up  to  l6  may  be  input. 

(F5.1) 

10 

Missile  Angles 
of  Attack 

XAL  Any  number  of  angles  of  attack  up  to  H8  may  be  input. 

(F5.1) 

11 

Missile  Angles 
of  Attack 

If  this  card  is  not  required,  leave  out  of  data  deck. 

12 

Missile  Angles 
of  Attack 

If  this  card  is  not  required,  leave  out  of  data  deck. 

TABLE  4 

OUTPUT  NOIffiNCLATUKE 


AL 

missile  angle  of  attack,  degrees 

CA 

total  axial  force  coefficient 

CDTOT 

total  missile  drag  coefficient 

CLALT 

of  the  tail 
a ' 

of  the  wing 

Li 

a 

of  wing  two 
a 

body  lift  coefficient 

CLALW 

CLALW2 

CLB 

CLI 

lift  loss  due  to  downwash 

CLTOT 

total  missile  lift  coefficient 

CLTT 

tail  lift  coefficient 

CLWT 

wing  lift  coefficient 

CM 

total  pitching  moment  coefficient 
gravity 

about  the  missile  center  of 

CN 

total  normal  force  coefficient 

CNB 

■  body  normal  force  coefficient  - 

CUT 

tail  normal  force  coefficient 

CNTD 

tail  normal  force  coefficient  due 
tion 

to  control 

surface 

deflec” 

CNW 

wing  normal  force  coefficient 

DELTA 

control  surface  deflection  angle, 

degrees 

VXM 

'  missile  flight  Mach  number 

XCPB 

body  center  of  pressure  location 
feet 

as  measured 

from  the 

nose , 

XCPT 

tail  center  of  pressure  location 
feet 

as  measured 

from  the 

nose , 

88 


TABLE  U 
(continued) 


XCPW 


ving  center  of  pressure  location  as  measured  from  the  nose, 
feet 

total  missile  center  of  pressure  location  as  measured  from 
the  nose,  feet 


XCP2 
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